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Preface 


This  report  describes  improvements  developed  for  foe  planetary  boundary 
layer  surface  wind  field  model  traditionally  used  by  the  U.S.  Army  Corps  of 
Engineers  for  hurricane  modeling.  Limitations  of  foe  model  are  also  described. 
The  upgraded  model  has  increased  flexibility  for  spatial  resolution  and  pressure 
profile  specification.  The  wind  fields  can  be  used  in  ocean  response  modeling, 
including  wave  and  surge  modeling  activities. 

This  study  was  authorized  by  Headquarters,  U.S.  Army  Corps  of  Engineers, 
under  foe  Coastal  Flooding  and  Storm  Protection  Area  of  foe  Coastal  Research 
Program,  Work  Unit  32683,  “Wind  Estimation  for  Coastal  Modeling.”  Techni¬ 
cal  Monitors  were  Messrs.  John  H.  Lockhart,  Jr.;  John  G.  Housley,  Barry  W. 
Holliday;  and  John  Saucier.  Ms.  Carolyn  M.  Holmes  of  foe  U.S.  Army  Engi¬ 
neer  Waterways  Experiment  Station  (WES),  Coastal  Engineering  Research 
Center  (CERQ,  was  foe  Program  Manager. 


The  study  was  conducted  under  Contract  No.  DACW39-93-C-0022  by 
Oceanweather,  Inc.  (OWI),  Cos  Cob,  Connecticut  The  report  was  prepared  by 
Dr.  Vincent  J.  Caidone  and  Messrs.  Andrew  T.  Cox  and  J.  Arthur  Greenwood, 
all  of  OWL  and  Dr.  Edward  F.  Thompson  of  foe  Coastal  Oceanography 
Branch  (COB),  Research  Division  (RD),  GERG  Dr.  Thompson  was  Principal 
Investigator  of  foe  research  work  unit  funding  this  study.  The  work  unit  was  Accession  Tor 
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At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


Availability  Codes 


Diet 


rvail  and/or 
Special 


The  conunts  of  Ms  report  an  mot  to  be  toed  for  adwnuing,  pubtication, 
or  promotional  pnrpam.  Citation  ef  trade  nomas  dots  not  conntntr  on 
(WWfaW  tmdonommi  tar  aonroMi  of  tko  uv  of  mdk  coototordol  nrnninrfi 


V 


□  □ 


Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


1  Multiply 

By . 

To  Obtain  I 

1  dogmas  (angto) 

0.01745329 

rattans  | 

|  knots  (international) 

0.5144444 

mows  par  sooond  | 

|  miss  (U.S.  nautical) 

1.852 

kiomotm  | 

vi 


1  Introduction 


Background 

The  unprecedented  destruction  of  parts  of  the  United  States  caused  by  hur¬ 
ricanes  Andrew  and  Iniki  last  summer  has  aroused  increased  interest  in  the 
wind  structure  of  tropical  cyclones  among  the  scientific  and  engineering  com¬ 
munities.  Unlike  most  past  destructive  storms,  much  of  the  loss  in  these  recent 
stomas  was  associated  with  direct  wind  damage.  While  simple  parametric 
tropical  cyclone  wind  models  remain  in  use  to  model  surface  winds  and  to 
provide  forcing  for  ocean  response  models,  a  few  numerical  vortex  boundary 
layer  models  based  upon  solution  of  foe  primitive  equations  of  motion  have 
emerged,  including  most  prominently  the  so-called  U.S.  Army  Corps  of  Engi¬ 
neers  (CE)  wind  model  (Cardone  et  al.  1992).  The  model  was  developed 
originally  at  New  York  University  in  the  early  1970’s  and  later  further  devel¬ 
oped  at  Oceanweather  Inc.  (OWI)  under  CE  support  in  1979.  Recently  the 
format  of  the  CE  model  was  modified  to  conform  with  die  requirements  of  the 
CE  Coastal  Modeling  System  (Thompson  1993). 

The  CE  numerical  model  was  reviewed  as  part  of  a  Workshop  on  Tropical 
Wind  Modeling  convened  at  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  on  24-23  March,  1992.  Invited  participants  in  the  review  were: 

Dr.  Wilson  A.  Schaffer,  Techniques  Development  Laboratory,  National 
Weather  Service  (NWS),  National  Oceanic  and  Atmospheric  Adminis¬ 
tration  (NOAA) 

Dr.  Mark  D.  Powell,  Hurricane  Research  Division  (HRD),  Atlantic  Oceano¬ 
graphic  and  Meteorological  Laboratory,  Environmental  Research  Labo¬ 
ratories,  NOAA 

Dr.  Muloit  B.  Mathur,  National  Meteorological  Center,  NWS,  NOAA 

Dr.  Vincent  J.  Cardone,  OWI 

The  workshop  stressed  the  relationship  between  surface  wind  modeling  and 
more  general  questions  of  foe  dynamical  and  thermodynamic  nature  of  tropical 
cyclones.  It  also  emphasized  the  need  to  carefully  evaluate  die  reliability  and 
representativeness  of  the  scant  surface  marine  wind  data  available  in  intense 
cyclones,  before  such  data  are  used  to  further  develop  and  validate  numerical 
models. 
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The  workshop  addressed  the  potential  for  further  development  of  existing 
numerical  models,  particularly  the  CE  model.  A  number  of  specific  research 
needs  for  improving  wind  models  were  identified  and  prioritized  >n  terms  of 
both  importance  to  ocean  response  modeling  and  feasibility  of  success.  These 
needs  are  summarized  more  fully  in  a  ‘‘white  paper*’  (Cardone  and  Thompson 
1992). 

Subsequent  to  the  workshop,  a  study  was  initiated  to  address  high  priority 
upgrades  to  the  CE  model  which  could  be  accomplished  with  the  limited  funds 
available.  The  following  tasks  were  chosen:  (1)  increase  resolution  and 
domain  of  the  nested  gad  system;  (2)  generalize  the  surface  pressure  specifi¬ 
cation.  The  enhancements  developed  for  the  CE  wind  model  are  the  subject  of 
this  report  The  enhancements  will  be  incorporated  into  the  Coastal  Modeling 
System  in  the  near  future. 


Previous  Studies 

The  CE  wind  model  has  been  used  mainly  to  provide  wind  fields  in  histori¬ 
cal  tropical  cyclones  to  drive  ocean  response  models  operated  in  a  hindcast 
mode  (surface  waves,  mixed  layer  currents,  storm  surge).  Those  wind  fields 
generally  provide  unbiased  ocean  predictions  when  used  to  drive  CE  and  OWI 
ocean  response  models  (e.g  Reece  and  Cardone  1982).  They  have  also  been 
used  to  drive  ocean  response  models  developed  by  other  scientists  indepen¬ 
dently,  wherein  CE  model  winds  have  also  repeatedly  been  shown  to  provide 
unbiased  hindcasts  (e.g.  Forristall  1980,  WAMDI  Group  1988,  Cooper  and 
Thompson  1989,  Ly  and  O’Connor  1991,  Grosskopf  et  al.  1991,  Mains  et  aL 
1992).  At  OWI  the  model  has  been  used  in  over  three  dozen  studies  to  drive 
ocean  response  models  to  establish  offshore  design  criteria  in  many  parts  of 
the  world  affected  by  tropical  cyclones. 

The  CE  model  has  been  extensively  used  for  both  ocean  wave  and  storm 
surge  modeling  for  CE  applications.  Abel  et  al.  (1989)  applied  the  model  to 
estimate  wave  statistics  due  to  hurricanes  in  the  Atlantic  Ocean  and  Gulf  of 
Mexico  during  1956-75.  Tracy  and  Hubertz  (1990)  estimated  waves  produced 
by  10  hurricanes  impacting  southern  California  during  1956-89.  Mark  and 
Scheffher  (1993)  describe  a  hurricane  surge  study  for  the  coast  of  Delaware. 

A  similar  approach  is  presently  being  applied  to  foe  entire  U.S.  Atlantic  Coast. 

The  generality  of  foe  CE  model  was  also  demonstrated  when  it  was  used  to 
provide  winds  to  test  the  third  generation  wave  model  (3GWAM)  (WAMDI 
Group  1988).  Winds  supplied  to  the  WAM  model  were  exactly  the  same  as 
winds  for  foe  subject  storms  (three  intense  Gulf  of  Mexico  hurricanes)  which 
had  been  used  in  previous  studies  to  drive  first  and  second  generation  models, 
and  which  had  been  used  by  other  investigators.  The  WAM  model  was  found 
to  provide  unbiased  and  skillful  wave  hindcasts  in  these  storms,  with  WAM 
using  its  own  calibration  of  source  terms  developed  completely  independently 
of  CE  winds.  The  same  tuning  on  WAM  has  also  been  shown  to  provide 
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neatly  perfect  hindcasts  in  severe  extratropical  storms  as  well  when  driven  by 
extremely  accurate  wind  fields  derived  by  direct  kinematic  analysis  of  wind 
measurements  (Cardone  et  al.  1994). 


In  many  of  the  studies  cited  above,  ocean  response  models  were  used  to 
evaluate  the  most  extreme  response  (storm  peak  winds,  waves,  surge  and  cur¬ 
rents)  in  a  storm  at  a  fixed  site.  In  general,  in  storms  in  which  die  assumed 
storm  pressure  profile  fits  the  actual  radial  distribution  well,  modeled  storm 
peaks  are  unbiased  in  the  mean  and  exhibit  scatter  index  of  IS  percent  or  less. 
The  method  is  less  successful  in  modeling  the  entire  spatial  and  temporal  dis¬ 
tribution  of  the  wind  field  in  such  storms.  There  are  some  storms  in  which 
even  the  storm  peaks  are  difficult  to  simulate,  where  the  storm  structure 
departs  from  die  simple  structure  implied  by  the  presumed  pressure  distribu¬ 
tion.  This  and  other  limitations  of  the  CE  model  are  described  in  more  detail 
in  Chapter  2. 

In  this  study,  two  limitations  of  the  CE  tropical  storm  wind  model  are 
addressed  and  remedied.  The  first  change  is  simply  die  addition  of  two  addi¬ 
tional  nests  to  die  grid  system  used  to  implement  the  numerical  vortex  model. 
This  change  provides  lower  truncation  errors  near  the  center  of  small  intense 
storms,  greater  resolution  near  die  vortex  center,  and  an  expanded  solution 
domain.  The  second  change  is  generalization  of  the  radial  surface  pressure 
profile  upon  which  the  surface  pressure  initialization  of  the  vortex  model  is 
based.  The  form  adopted  also  allows  die  specification  of  profiles  with  two 
maxima  in  the  radial  pressure  gradient.  These  changes  are  described  in 
Chapter  3.  A  summary  is  given  in  Chapter  4. 
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2  Existing  Model  Limitations 


The  CE  wind  model  developed  by  Caidone  et  al.  (1992)  has  proved  to  be  a 
powerful  tool  in  ocean  response  modeling.  However,  the  model,  developed  in 
1979,  includes  a  number  of  limitations.  In  light  of  enhanced  computing  power 
now  available  and  the  increasing  field  measurements  and  understanding  of 
tropical  storm  behavior,  it  is  timely  to  review  tire  model  limitations.  Limita¬ 
tions  of  the  CE  model  may  be  described  in  three  basic  categories:  physics, 
initialization,  and  numerics. 


Physics 

The  CE  model  evolved  from  the  model  of  Chow  (1971)  who  solved  tire 
momentum  equations  of  an  integrated  boundary  layer  flow  for  a  boundary 
layer  of  constant  depth.  The  vertical  friction  force  was  taken  parallel  to  the 
wind  relative  to  foe  earth.  Horizontal  friction  was  also  considered.  The  equa¬ 
tions,  however,  were  solved  numerically  on  a  nested  cartesian  grid  system 
centered  on  the  vortex  and  translating  at  constant  velocity  with  the  vortex. 

The  steady  solution  in  foe  moving  coordinate  system  referred  back  to  foe  earth 
yielded  qualitatively  realistic  boundary  layer  wind  patterns.  The  solution 
included  supergradient  flows  inside  foe  radius  of  maximum  gradient  wind  and 
a  decrease  in  the  radius  of  maximum  wind.  It  also  included  an  asymmetric 
wind  distribution  with  maxima  in  the  right  front  quadrant  for  a  typical  super¬ 
position  of  a  symmetric  vortex  and  ambient  gradient,  and  a  boundary  layer 
convergence  pattern  consistent  with  observed  patterns  of  convection  in  typical 
storms. 

Shapiro  (1983)  solved  the  same  slab  momentum  equation  as  Chow  (1971) 
but  used  a  truncated  spectral  analysis  in  cylindrical  coordinates,  in  order  to 
allow  a  more  convenient  separation  of  the  role  of  linear  and  non-linear  asym¬ 
metric  effects  in  foe  boundary  layer  flow.  Chow’s  model  and  solution  method 
provide  foe  same  patterns  as  that  of  Shapiro’s  model  except  that  inside  the 
radius  of  maximum  wind  truncation  errors  are  larger  than  for  the  spectral  solu¬ 
tion.  As  a  consequence,  Chow’s  model  may  slightly  overestimate  foe  degree 
of  supergradient  flow  inside  the  eye.  These  studies  straw  that  foe  essential 
physics  governing  the  boundary  layer  flow  are  included  in  Chow’s  and 
Shapiro’s  models.  The  main  physical  processes  missing  are  foe  feedback  of 
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the  convection  (induced  in  part  by  the  modelled  convergence  field)  on  the 
wind  field,  and  strong  non-steady  effects  (for  example  rapid  deepening  of  the 
vortex  in  the  moving  frame)  which  may  cause  even  the  overlying  vortex  to  be 
unbalanced. 

The  CE  model  is  derived  directly  from  Chow’s  formulation  and  ures 
Chow’s  numerical  solution.  Several  improvements  to  Chow’s  solution  were 
made  to  insure  that  it  not  only  gives  qualitatively  realistic  wind  fields,  but  also 
provides  a  quantitatively  correct  surface  stress  vector  distribution,  from  which 
the  model  diagnoses  winds  within  die  surface  boundary  layer  as  well.  The 
main  enhancements  to  Chow’s  model  are  in  the  inclusion  of  a  similarity 
boundary  layer  formulation  relating  vertically  integrated  flow  to  the  surface 
drag  (magnitude  and  direction),  adoption  of  more  realistic  boundary  layer 
depths  than  considered  by  Chow  or  Shapiro,  consideration  of  the  effects  of 
boundary  layer  stratificaL  and  variable  surface  roughness  (expressed  in  terms 
of  wind  alone  with  no  sea  state  effects  considered),  and  incorporation  of 
greater  flexibility  in  die  specification  of  the  imposed  pressure  distribution  of 
the  vortex  over  the  possibilities  considered  by  Chow. 

The  CE  model  was  dc. eloped  with  a  secondary  objective  to  provide  winds 
over  inland  lake  surfaces  and  land  surfaces  of  arbitrary  roughness.  The  theo¬ 
retical  development  of  this  part  of  the  model  met  with  less  success  than  the 
over-water  treatment.  A  simplified  equilibrium  boundary  layer  approach  was 
adopted  which  ignores  the  adjustment  of  the  planetary  boundary  layer  (PBL) 
wind  field  across  discontinuities  of  roughness.  Thus,  while  die  model  valida¬ 
tion  against  winds  measured  over  land  indicated  good  agreement  when  the 
wind  fetch  was  over  a  homogeneous  roughness,  little  is  known  about  the  effect 
on  ocean  response  modeling  associated  with  failure  in  the  CE  model 
(probably)  to  resolve  small  scale  PBL  wind  changes  downwind  of  abrupt 
changes  in  roughness  (e.g.  the  coast). 


Initialization 

The  model  is  generally  applied  with  boundary  layer  height  in  the  range  of 
500  m  -  650  m,  slightly  unstable  stratification,  a  Chamock  type  surface  rough¬ 
ness  formulation  (Chamock  constant  0.035  with  Karman  constant  0.35),  and  a 
value  of  unity  for  the  Ekman  scale  height  parameter.  This  combination  pro¬ 
duces  unbiased  surface  winds  over  the  open  sea  when  the  model  is  applied  to 
real  storms  and  validated  against  measured  surface  wind  time  histories 
obtained  by  calibrated  instruments  (e.g.  NOAA  buoys,  offshore  rigs). 

The  pressure  field  is  generally  described  as  the  superposition  of  the  pressure 
gradient  computed  from  the  exponential  pressure  profile  form  for  the  symmet¬ 
ric  part  of  the  vortex: 
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p(r )  ■  P,  *  <P--Pje  ' 


(1) 


where 

p(r)  =  pressure 

p„  =  central  pressure  (at  the  eye) 

pm  =  axisymmetric  ambient  pressure  (far  field  pressure) 

=  scaling  radius 
r  »  radius 

and  an  uniform  ambient  gradient  given  by 

fix  TT  -  -  I  (2) 


where 


/«  Coriolis  parameter 
X  =  unit  vector  in  die  vertical  direction 
Vg  =  ambient  uniform  geostrophic  flow 
p  =  mean  air  density 
Vpa  s  uniform  ambient  pressure  gradient 


This  pressure  initialization  scheme  (it  is  also  a  boundary  condition  since  the 
model  is  solved  to  a  steady  state  solution)  often  provides  a  very  realistic  simu¬ 
lation  of  the  actual  pressure  field  about  a  tropical  cyclone.  However,  in  some 
storms  the  actual  pressure  field  departs  from  this  simple  picture  in  several 
possible  ways.  Often,  particularly  as  a  tropical  cyclone  enters  die  mid¬ 
latitudes,  the  ambient  pressure  field  is  inhomogeneous.  The  effect  is  especially 
evident  if  die  tropical  cyclone  begins  to  interact  with  a  frontal  system  or  an 
extratropical  cyclone  or  both.  Within  the  tropics,  some  storms  have  been 
shown  (Holland  1980)  to  follow  the  more  general  form: 
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p(r)  *  Pe 


(P.  -  Pj  * 


where 

B  =  constant  in  the  general  range  05-25 

Finally,  in  some  storms  the  radial  pressure  profile  in  the  inner  core  is  more 
irregular  than  either  of  the  above  forms,  with  a  shape  which  implies  two 
maxima  in  the  radial  pressure  gradient,  accompanied  by  two  distinct  maxima 
in  the  wind  speed.  Willoughby  (1990)  and  Black  and  Willoughby  (1992)  have 
described  the  tendency  for  “concentric  rings”  in  the  radial  wind  distribution  to 
be  a  fairly  typical  characteristic  of  intense  tropical  cyclones.  The  rings  appear 
to  be  related  intimately  to  storm  intensity  evolution.  For  example,  Figure  1 
shows  the  evolution  of  concentric  rings  in  Hurricane  Gilbert  (1988)  over  a 
six-day  period.  In  this  storm,  foe  CE  model  might  be  expected  to  provide 
reasonably  accurate  wind  fields  in  foe  initial  stage  of  vortex  development  and 
intensification  between  September  1 1-13,  but  it  would  fail  to  model  foe  com¬ 
plicated  double  maxima  structures  later.  The  impact  mi  ocean  response  model¬ 
ing  of  this  failure  to  model  concentric  rings  is  unknown. 


Numerics 


The  CE  model  is  computationally  demanding.  For  example,  computational 
considerations  drove  foe  decision  to  presolve  foe  boundary  layer  model  for 
spatially  homogeneous  (constant)  boundary  layer  height  and  stability  and  to 
use  a  table  look-up  procedure  for  foe  drag  coefficient  during  foe  matching  of 
the  solution  toward  steady  state.  To  relax  the  constraint  of  constant  boundary 
layer  height  and  stability  would  greatly  increase  computer  time,  unless  a  more 
efficient  integration  scheme  could  be  found.  The  minimum  grid  spacing  on 
the  inner  nest  of  the  solution  grid  (as  opposed  to  the  target  grid)  is  5  km, 
which  is  a  bit  too  large  to  resolve  details  of  the  wind  field  near  foe  center  in 
very  tight  storms.  (For  example,  as  Hurricane  Andrew  approached  the  smith 
Florida  coast,  foe  radius  of  maximum  wind  was  only  about  11  km.)  The  grid 
spacing  is  also  not  sufficient  to  resolve  boundary  layer  adjustments  near  rough¬ 
ness  discontinuities,  though  that  physical  process  is  not  presently  incorporated 
in  the  numerical  model.  Further,  the  grid  spacing  is  too  coarse  to  resolve 
detailed  gradients  of  wind  over  inland  bays  and  estuaries.  Limitations  in  tem¬ 
poral  resolution  are  less  serious,  within  the  constraint  of  the  steady-state 
model,  since  the  “time  step”  of  the  windfields  is  simply  the  temporal  resolu¬ 
tion  of  the  storm  track.  That  temporal  resolution  can  be  refined  within  reason 
(say  to  intervals  of  13  minutes  or  so)  without  significant  computational  cost 
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TANGENTIAL  NINO  (•/«! 


Figure  1 .  Temporal  changes  in  the  azimuthally  averaged  wind  for  NOAA  reconnaissance 
flights  into  Hurricane  Gilbert;  changes  are  normalized  to  a  6-hr  time  interval  (from 
Black  and  Willoughby  (1992)) 
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Summary  of  Limitations 


The  main  model  limitations  described  in  each  category  above  may  be  listed 
as  follows: 

a.  Physics. 

(1)  Decoupling  of  boundary  layer  from  full  vortex  dynamics,  pre¬ 
cluding  mutual  adjustment  of  pressure  and  wind  fields;  and  feed¬ 
back  of  convective  scale  effects  on  wind  field. 

(2)  Simplified  PBL  theory  (e.g.  constant  Ekman  scale  height). 

(3)  Extrapolation  of  Chamock  roughness  to  extreme  wind  speeds,  with 
no  sea  state  dependence. 

(4)  No  boundary  layer  adjustment  across  roughness  discontinuity. 

b.  Initialization. 

(1)  Constant  and  homogeneous  boundary  layer  height 

(2)  Constant  and  homogeneous  stratification. 

(3)  Relatively  simple  pressure  specification: 

(a)  Exponential  pressure  profile  provides  only  one  radius  of 
maximum  wind  (no  concentric  rings). 

(b)  Pressure  profile  may  be  inadequate  even  for  unimodal 
maximum  gradient  pressure  distributions. 

(c)  Homogeneous  ambient  linear  pressure  gradient. 

c.  Numerics. 

(1)  Practical  limit  of  spatial  resolution  to  5  km  may  be  inadequate  for 
very  tight  storms. 

(2)  Large  number  of  iterations  (800)  required  for  each  steady  state 
configuration,  or  snapshot  (using  the  terminology  of  Cardone  et  al. 
1992). 

In  this  study,  two  of  die  above  limitations  are  addressed  and  remedied,  as 
described  in  the  next  section. 
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Increased  Resolution 

The  CE  program  consists  basically  of  two  main  programs  (Cardone  et  al. 
1992).  The  first,  SNAP,  solves  the  numerical  vortex  model  on  a  nested  grid  a 
number  of  times  to  represent  die  storm  wind  field  at  discrete  times  within  an 
event,  thereby  producing  a  number  of  snapshot  wind  fields  on  a  nested  grid. 
SNAP  also  writes  the  snapshots  to  a  file  for  use  by  the  second  main  program 
HIST,  which  among  other  functions,  interpolates  the  nested  grid  solutions  to 
hourly  intervals  and  then  interpolates  the  winds  to  an  output  or  target  grid 
(typically  that  of  an  ocean  response  model). 

The  nested  grid  consists  of  five  square  21  by  21  grid  point  arrays.  The  grid 
spacing  increases  by  a  factor  of  two  from  nest  to  nest.  In  the  program  input, 
the  user  specifies  a  desired  spacing  of  die  inner  nest  (die  variable  DX  in  name- 
list  NAME3  of  SNAP).  For  the  default  value  of  5  km,  die  grid  spacing  in  the 
coarsest  nest  becomes  80  km  and  die  attire  grid  covers  an  area  of  (1,600  km1). 
While  die  existing  CE  program  allows  users  to  set  DX  smaller  than  5  km,  this 
is  not  recommended  since  the  grid  coverage  shrinks  commensurately.  The 
simplified  boundary  condition  applied  on  die  outer  boundary  of  die  outermost 
nest  becomes  increasingly  tenuous  as  die  gridded  domain  shrinks. 

The  upgraded  program  allows  the  use  of  up  to  seven  nests.  However,  the 
user  may  specify  the  number  of  nests  (from  three  to  seven)  in  a  given  run. 

The  new  parameter  INSIDE  is  used  to  specify  die  number  of  active  nests.  It 
designates  which  is  the  finest  active  nest,  where  nests  are  numbered  from  1  to 
7  going  from  finest  to  coarsest  For  example,  INSIDE  =  1  activates  all  seven 
nests,  and  INSIDE  =  2  activates  only  nests  2  through  7.  The  grid  spacing  of 
nest  1,  the  innermost  nest  is  specified  as  before  with  variable  DX,  regardless 
of  whether  or  not  the  nest  is  active.  The  relationships  between  INSIDE,  active 
nests,  and  spacing  of  die  finest  active  nest  are  summarized  in  Table  1.  The 
default  value  of  DX  is  2  km.  If  all  seven  nests  are  exercised,  the  execution 
time  per  snapshot  is  roughly  four  times  as  long  as  die  existing  CE  model.  For 
this  case  the  number  of  iterations  on  the  inner  nest  is  set  to  the  default  value 
of  3200. 
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Table  1 

Effect  of  Neat  Activation  Parameter,  INSIDE 

INSIDE 

Active  NaH 

Spacing  of  Hnoat  Active  Neat  | 

1 

1  -7 

OX 

2 

2-7 

2  DX 

3 

3-7 

4  OX 

4 

4-7 

8  DX 

5 

5-7 

16  DX 

Hie  program  with  die  new  nesting  was  tested  in  two  ways  using  the  Hurri¬ 
cane  Camille  snapshots  as  test  cases.  First,  vortex  model  winds  were  produced 
by  the  GE  model  (and  its  OWI  equivalent)  for  the  case  of  DX  =  8  km.  Then, 
the  same  SNAP  inputs  were  used  to  generate  winds  with  the  new  code  for  the 
case  DX  =  2  km,  INSIDE  =  3,  which  provides  die  equivalent  number  of  nests 
and  inner  nest  grid  spacing  as  the  CE  model  run.  Winds  produced  by  the  two 
alternative  programs  were  interpolated  to  a  target  grid  covering  die  Gulf  of 
Mexico  (nominal  spacing  of  0.2  deg1),  compared  and  found  to  agree  to  within 
roundoff  error  of  the  VAX  computer  used  for  these  tests. 

The  second  test  compared  winds  for  Camille  produced  by  die  new  code  for 
the  case  DX  =  2  km  and  INSIDE  =  1;  that  is,  all  seven  nests  are  live  with 
inner  nest  grid  spacing  of  2  km,  with  winds  produced  by  die  CE  program  with 
DX  =  5  km.  These  results  are  shown  in  Appendix  A  which  gives,  at 
30-minute  intervals,  die  maximum  scalar  wind  speed,  the  location  and  die  cor¬ 
responding  wind  speeds  and  directions,  and  die  same  data  for  die  maximum 
vector  wind  difference  magnitude.  The  results  (see  also  Figure  2)  indicate  that 
the  largest  differences  (scalar  differences  of  up  to  about  7  m/s),  occur  inside 
the  eyewall,  where  truncation  errors  on  the  5-km  solution  are  expected  to  be 
large  for  an  intense  tight  vortex  such  as  Camille.  Maximum  scalar  wind  speed 
differences  in  the  area  of  the  eyewall  are  generally  less  than  1  m/sec.  How¬ 
ever,  maximum  vector  difference  magnitudes  of  up  to  9  m/sec  were  observed 
occasionally  in  the  vicinity  of  the  eyewall  reflecting  a  tendency  for  the  wind 
direction  on  die  2-km  solution  to  be  turned  systematically  in  the  direction  of 
less  inflow,  by  up  to  10  deg  from  the  5-km  solution.  The  2-km  solutions  are 
no  doubt  the  more  accurate  solutions. 


1  A  table  of  factors  for  converting  non -SI  units  of  measurement  to  SI  units  is  presented  on 
page  vi. 
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Figure  2.  Distribution  of  maximum  wind  speed  differences  between  5-nest 
and  7-nest  model  runs  for  Hurricane  Camille 


Generalized  Pressure  Specification 

Pressure  profile  form 

The  upgrade  to  die  pressure  specification  uses  a  generalized  form  of 
Holland’s  (1980)  exponential  pressure  profile 


-  -/M* 

Pir)  -  p0  *  Y.  4P,«'r/ 


where 


n  =  number  of  components 

dpi  =  pressure  anomaly  for  the  i'th  component  (dpt+dp2+...+dp,,=p„-pJ 
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Rp  m  scaling  radius  for  the  i’th  component 
Bj  =  Holland's  B  coefficient  for  the  i’th  component 
The  corresponding  tangential  and  radial  pressure  gradients  are: 


dp 
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(i tangential ) 
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ft 
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(5) 


(radial) 


The  CE  model  was  modified  to  provide  options  for  a  single  or  double  expo¬ 
nential  component  (1*  «  7  or  n  *  2).  This  form  allows  the  specification  of 
pressure  profiles  with  two  separate  maxima  in  the  radial  pressure  gradient, 
though  the  mere  form  does  not  guarantee  two  maxima.  For  example,  the  sum 
of  two  exponentials  also  allows  the  modeling  of  pressure  profiles  which  have 
only  one  maximum  but  with  shapes  very  different  from  those  predicted  by  the 
single  exponential  even  with  the  variable  B  included. 

Incorporation  of  this  model  into  the  seven-nest  version  of  the  program  led 
to  extensive  changes  to  program  SNAP,  in  particular  to  Subroutines  GRAD 
and  PXYM  and  further  changes  to  namelist  NAME3,  as  documented  in 
Appendix  B.  One  immediate  consequence  of  this  model  is  that  the  quadrantal 
specification  of  profile  parameters  allowed  in  tire  single  exponential  form  with 
B  ss  7  is  lost  To  retain  this  option  in  the  CE  model  portfolio,  two  versions  of 
the  upgraded  program  were  developed  as  follows: 

SNAP_ADC.7NE  and  HIST_ADC.7NE  -  This  version  only  upgrades 
the  current  CE  program  (which  includes  quadrantal  variation  of  parame¬ 
ters  for  a  single  exponential  with  B  =  7)  to  incorporate  the  additional 
nests. 


SNAP_HOL.7NE  -  This  version  upgrades  tire  current  CE  model  to 
allow  both  7  nests  and  tire  generalized  pressure  specification  scheme, 
but  without  quadrantal  variation.  Note,  however,  that  asymmetry  in  the 
pressure  field  is  still  modeled  through  superposition  of  the  vortex  pres¬ 
sure  field  and  the  background  steering  gradient  The  background  pres¬ 
sure  gradient  is  required  to  be  homogeneous  (that  is,  the  parameter 
ST12  is  eliminated). 


Modified  outflow 

In  early  tests  of  the  upgraded  CE  program  with  hypothetical  snapshot  inputs 
estimated  roughly  to  apply  to  several  stages  of  Hurricane  Gilbert  (inputs  for 
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Gilbert  snapshots  were  derived  more  rigorously  as  described  below),  it  was 
demonstrated  that  the  program  could  produce  the  pattern  of  annular  concentric 
wind  maxima.  However,  it  was  noticed  that  the  inflow  characteristics  of  the 
model  appeared  to  have  changed  somewhat  from  the  unimodai  mode!.  In  die 
standard  CE  model,  Subroutine  OUTFLOW  serves  to  remove  8  degrees  of 
inflow  from  the  snapshot  solution  throughout  the  domain  to  compensate 
approximately  for  inflow  believed  to  be  spuriously  introduced  through  the 
numerical  solution.  This  spurious  inflow  is  revealed  by  solving  the  model  for 
a  motionless  vortex  with  all  friction  terms  deactivated  and  comparing  the 
modeled  wind  direction  to  the  purely  circular  flow  expected  of  a  vortex  in 
gradient  balance. 

The  frictionless,  motionless  vortex  test  was  repeated  with  the  upgraded 
model  for  the  series  of  nine  snapshots  indicated  in  Table  2.  For  the  unimodai 
case  (Case  1)  and  the  bimodal  cases  with  B  -  1  for  both  exponentials  (cases  2 
and  6),  the  maximum  inflow  (which  tends  to  occur  just  outside  the  wind 
maxima)  averages  75  deg.  As  the  value  of  B  increases,  however,  there  is  a 
proportional  increase  in  tile  inflow.  In  the  bimodal  profile  with  B  varying 
between  the  two  exponentials,  the  inflow  varies  with  radius.  Subroutine  OUT¬ 
FLOW  was  modified  to  compensate  for  dependence  of  spurious  inflow  on  B 
(Table  3). 


I  Table  2 
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Empirical  Correction  of  Inflow  Angle 

tjpfirundal 

Reduction  Applied  To  Mow  Anglo  | 

UnimodaJ  0-  1 

. 

Unimode)  Sul 

80  ! 

^ - 

nmooii 

•M  ] 

Specification  of  pressure  parameters 

Single  Exponential  Profile:  B  =  1.  In  the  standard  CE  model,  the  pressure 
profile  may  be  specified  in  basically  two  ways.  The  most  fundamental  way  is 
to  fit  die  profile  to  sea  level  pressure  measurements  available  at  different  radii 
at  a  given  time,  or  transformed  from  time  to  space  using  single  station  data 
acquired  at  a  station  in  die  path  of  die  storm.  There  are  several  examples  of 
this  approach  as  applied  to  historical  U.S.  Gulf  of  Mexico  and  East  Coast  hur¬ 
ricanes  in  Graham  and  Hudson  (1960).  The  eye  pressure  may  be  prescribed,  if 
it  is  known,  for  a  more  accurate  fit,  or  the  eye  pressure  may  be  extrapolated 
from  a  fit  determined  exclusively  from  data  outside  the  center.  A  simplifica¬ 
tion  to  this  procedure  is  often  followed  for  oceanic  storms,  for  which  eye  pres¬ 
sure  may  be  known  from  aircraft  dropsonde  data,  far  field  pressure  is  estimated 
from  weather  maps,  and  a  few  estimates  of  pressure  at  various  radii  about  the 
storm  are  known  from  ship  or  island  station  synoptic  reports.  Then,  die 
unknown  parameter  scale  radius  may  be  estimated  from  Equation  1  as  follows 
for  each  such  report  and  an  average  or  weighted  average  of  the  estimates  taken 
to  represent  the  storm  profile  at  map  time: 


-  r  In 


Pto  -  P0 


VP--P., 


(6) 


If  there  are  insufficient  pressure  data  but  eye  pressure  is  known  and  an 
estimate  of  the  radius  of  maximum  wind,  Rm,  is  known  (e.g.  from  aircraft  vor¬ 
tex  message  reports  filed  upon  penetration  of  the  eye  and  assuming  that  Rm  at 
flight  level  is  the  same  as  A,  at  the  surface,  or  more  crudely  from  radar  or 
satellite  eye  diameter  estimates),  Rp  may  be  estimated  directly  from  Rm  using 
the  average  relationship  found  between  these  two  variables  by  the  vortex 
model. 
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Single  Exponential  Profile:  Variable  B.  Variants  of  these  same  two 
approaches  may  be  followed  to  estimate  the  parameters  of  the  generalized 
iininwviai  model,  for  which  the  additional  parameter  B  must  also  be  specified. 
Again,  if  there  are  sufficient  pressure  data,  foe  profile  may  be  fitted  directly. 
For  example,  Figure  3  shows  foe  screen  display  of  a  PC-based  interactive 
system  developed  at  OWI  using  a  commercially  available  (dotting/statistical 
analysis  software  package.  The  pressure  data  are  composited  (see  window  in 
upper  right  hand  comer  of  screen)  as  a  function  of  radius  in  a  Smith  China  Sea 
typhoon  from  reduced  (from  flight  level)  aircraft  pressures  and  pressures 
reported  by  ships  within  a  3 -hr  time  window  of  analysis  time.  The  aircraft 
also  provided  estimates  of  eye  pressure  (note  the  two  conflicting  estimates  at 
the  origin  in  the  lower  window  of  Figure  3).  Far  field  pressure  was  estimated 
from  weather  maps.  The  best  fit  shown  is  for  pc  =  968  mb,  B  =  0.8,  and  Rp  = 
(14  nm)**.  The  window  at  the  upper  left  of  foe  screen  compares  the  azimuth- 
ally  averaged  solution  for  foe  model  surface  wind  (downloaded  from  a  run 
made  on  a  VAX)  and  reduced  aircraft  and  ship  reports  of  wind. 

The  second  approach  is  more  indirect  and  emphasizes  the  use  of  aircraft 
wind  data.  In  recent  years  such  data  have  become  quite  accurate  after  foe 
introduction  of  inertial  navigation  systems,  onboard  processing  and  the  avail¬ 
ability  of  coded  messages  (so-called  supplementary  or  peripheral  flight  level 
winds)  containing  measured  winds  at  flight  level  outside  foe  eye.  Figure  1 
shows  foe  complete  analysis  carried  out  by  Black  and  Willoughby  (1992)  of 
flight  data  acquired  over  foe  main  lifetime  of  Hurricane  Gilbert  in  the  Carib¬ 
bean  Sea  and  Gulf  of  Mexico.  These  curves  show  12  separate  radial  profiles 
of  the  azimuthally  averaged  flight  level  (700  mb  or  830  mb)  tangential  wind 
speed  composited  from  flight  legs  near  foe  indicated  times.  Most  of  the  wind 
profiles  exhibit  two  distinct  wind  maxima.  On  foe  assumption  that  the  azi¬ 
muthally  averaged  flow  is  approximately  in  gradient  balance  with  the  axisym- 
metric  pressure  field,  the  pressure  profile  associated  with  double  concentric 
wind  maxima  should  also  exhibit  two  local  maxima  in  the  pressure  gradient 
Therefore  they  might  be  fitted  by  a  double  exponential  profile.  Even  those 
profiles  which  do  not  exhibit  two  distinct  peaks,  such  as  those  in  panels  b,  d, 
and  e,  exhibit  atypical  shapes  for  tropical  cyclones,  with  a  single  maximum 
and  broad  regions  with  little  or  no  change  of  wind  speed  with  radius.  Never¬ 
theless,  we  have  selected  five  of  these  cases  to  illustrate  foe  fitting  of  foe 
single  exponential  profile  using  aircraft  wind  data  and  eye  pressure.  Parame¬ 
ters  are  defined  in  Table  4  and  results  from  the  fitting  process  are  given  in 
Table  5. 

The  fitting  method  is  basically  a  systematic  search  of  many  possible  solu¬ 
tions  of  foe  single  exponential  for  that  solution  whose  radial  distribution  of 
implied  gradient  wind  provides  a  close  match  to  the  location  and  magnitude  of 
the  azimuthal  average  flight  level  wind.  This  searching  program  (implemented 
in  a  preprocessing  program  called  lEYEWALL.nOL)  requires  foe  input  infor¬ 
mation  listed  in  Table  3.  Additional  documentation  is  given  in  Appendix  B. 
The  searching  program  fixes  foe  profile  anomaly  parameter  (p„  -  pj,  loops 
through  possible  values  of  foe  scale  radius  Rp  (from  Rm  to  2RJ  and  B  (from 
0.3  to  2.32)  and  finds  the  pressure  profile  whose  gradient  wind  simultaneously 
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Oceanweather  Tropical  System  Analysis 
Surface  Winds  and  Pressures  Estimated  from 
Reccon.  Vortex  and  Periperal  Data  Messages 


Typhoon  WAYNE86 
86090212  +/-  3hrs 


0  60  120  180  240 

Radial  Distance  from  Center  (Nml) 


Figure  3.  Example  of  Oceanweather  tropical  storm  analysis 
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Table  4 

Parameter  Definition*  for  Fitting  Single  Exponential  Profile 


Obearved  rectus  of  marimum  wind 


ASn*iWle^f  IWeQaO  BfijafllW  M^n  MVw  WflQ  H  n. 


AzimuttaNy  nwwgid  tangential  light  level  wind  at  150  km  from 


Seals  radhn  of  tingle  exponendaf  proflie 


Sham  preeetaa  anomaly  parameter  of  tingle  anponamial  proMe 


HotoKfs  * 


Hadkja  of  maximum  nnaient  of  Mad  mature  oro*a 


Maadmum  gradent  wind  of  Mad  prasturs  proSa 


Gradent  wind  of  Mad  protie  150  km  from  cantor 


Table  5 

Generalized  Single  Exponential  Profile  Fits  to  Selected  Hurricane  Gilbert 


Output 


&  Sr 


30 

20 

107 

1.50 

67.2 

30 

16 

124 

1.41 

70.3 

36 

97 

61 

1.00 

39.3 

38 

92 

62 

1.12 

42.2 

37 

85 

61 

1.12 

41.9 

best  matches  the  observed  wind  speeds  at  Rm  and  at  150  km  in  terms  of  abso¬ 
lute  difference.  For  example,  in  Case  3  of  Table  5  die  selected  profile  gradient 
wind  (not  shown)  is  within  1  m/s  of  the  Vm  of  65  m/s  at  Rm.  However  the 
searching  program  {daces  die  absolute  profile  maximum  of  67.2  m/s  at  a  radius 
of  20  km,  and  fails  to  maintain  the  observed  broad  region  of  little  change  in 
wind  speed  between  50  and  150  km,  resulting  in  a  profile  wind  speed  of  only 
21  m/s  at  150  km,  about  10  m/s  lower  than  measured.  Cases  7,  8  and  9  are 
somewhat  more  successful.  The  parameter  B  varies  between  1.0  and  1 J  for 
these  fits,  or  in  the  same  general  range  reported  by  Holland  for  a  single 
exponential. 

Double  Exponential  Profile.  While  it  is  conceivable  that  there  may  be 
sufficientiy  voluminous  and  accurate  pressure  data  in  some  tropical  cyclones  to 


18 


Chapter  3  Upgraded  Modal 


attempt  to  directly  fit  a  double  exponential  profile  to  measurements  of  surface 
pressure  as  a  function  of  radius,  we  have  not  attempted  to  construct  such  a 
data  set  and  perform  such  a  fit.  We  did  try,  without  success,  to  develop  useful 
fits  to  the  profile  from  just  the  total  storm  anomaly,  estimates  of  the  two  radii 
of  maximum  wind  and  a  single  pressure  along  the  profile  in  the  region 
between  the  two  maxima.  However,  a  generalization  of  the  searching  algo¬ 
rithm  described  above  for  a  single  exponential  has  met  with  some  success. 

The  searching  algorithm  (called  2EYEWALL.HOL)  as  applied  to  a  double 
exponential  is  documented  in  Appendix  B.  Parameters  involved  are  defined  in 
Table  6  and  Figure  4.  The  searching  program  fixes  the  total  storm  anomaly, 
assumes  the  scale  radius  for  the  inner  exponential  is  equal  to  the  observed 
radius  of  maximum  wind  of  the  inner  maximum,  or  ring,  and  loops  through 
ranges  of  the  outer  scale  radius,  R#  (from  R^  to  2Rm 2),  inner  and  outer  B 
(from  0.5  to  2.52)  and  the  ratio  of  dpjdp2  (from  1/8  to  8).  The  algorithm 
seeks  the  combination  whose  pressure  profile  provides  a  gradient  wind  profile 
that  has  maxima  at  R*,  and  R^  within  ±1  m/s  of  observed  and  which  maxi¬ 
mizes  the  following: 

VtJ  *  ^  -  2^  (7) 


where 

=  gradient  wind  at  (R^+R^ft 

If  it  succeeds  in  finding  such  a  profile  it  checks  that  the  wind  at  150  km  is 
lower  than  the  outer  maximum,  V^,  and  if  it  is,  prints  the  solution.  If  these 
conditions  are  not  met,  another  cycle  is  attempted.  The  matching  criterion  is 
relaxed  to  ±2  m/s,  this  time  requiring  that  foe  wind  at  (R^+R^IZ  is  less  than 
the  winds  prescribed  at  R*,  and  R^,  and  minimizing  the  wind  at  150  km.  The 
program  also  prints  foe  profile  parameters  for  the  selected  profile.  If,  after  the 
second  cycle  foe  program  still  does  not  find  a  successful  fit,  it  prints  the 
closest  fit  found  in  each  cycle. 

Table  7  shows  the  results  of  the  application  of  this  searching  algorithm  to 
all  12  of  foe  azimuthal  average  tangential  flight  level  wind  profiles  in  Hurri¬ 
cane  Gilbert  derived  by  Black  and  Willoughby  (1992).  The  double  exponen¬ 
tial  appears  to  require  large  values  of  B  to  resolve  two  distinct  peaks  in  the 
radial  profile  of  gradient  wind,  at  least  for  most  of  these  cases.  Table  8  com¬ 
pares  the  location  and  magnitude  of  the  double  wind  maxima  derived  from  the 
fitted  profile  to  those  observed.  Figure  5  compares  the  fitted  and  observed 
radial  profiles  of  pseudo-gradient  wind.  The  inner  ring  is  usually  fitted  very 
closely.  In  1 1  out  of  12  of  the  cases  a  distinct  outer  wind  maximum  is 
resolved,  and  it  is  usually  placed  within  ±20  km  of  the  observed  maximum.  In 
9  out  of  those  11  cases,  tire  maximum  is  within  about  ±2  m/s  of  that  observed. 
Case  8  is  the  poorest  fit,  but  in  practice  Case  7  and  Case  8  are  so  close  in  time 
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Tables 

Parameter  Definitions  for  Fitting  Double  Exponential  Profile 

Perea** 

- y 

Definition 

Radius  of  maxim  urn  wind,  inner  ring 

FL, 

Radius  of  maximum  wind,  outer  ring 

v„, 

Azimuthal iy  averaged  tangential  flight  level  wind,  inner  ring 

Azimuthal ly  averaged  tangential  flight  level  wind,  outer  ring 

a. 

Scale  radius,  inner  exponential 

Scale  radius,  outer  exponential 

dp, 

Pleasure  anomaly,  inner  exponential 

dp. 

Pressure  anomaly,  outer  exponential 

B, 

Hotantfs  B,  inner  exponential 

B, 

Hotoxfs  B,  outer  exponential 

R* 

Radius  of  maximum  gradient  wind  of  fitted  profile,  inner  ring 

R* 

Radius  of  maximum  gradient  wind  of  fitted  profile,  outer  ring 

V.' 

Maximum  gradnnt  wind  of  fitted  profile,  inner  ring 

V* 

Maximum  gradtont  wind  of  fitted  profile,  outer  ring 

(in  fact  they  are  derived  from  the  same  flight)  that  Case  7  may  be  used  to  rep¬ 
resent  this  phase  of  die  storm  history. 


Sample  Runs 

The  upgraded  program,  including  seven  nests  and  die  generalized  pressure 
specification,  has  been  applied  to  provide  sample  wind  fields  on  target  grids 
using  as  test  input  the  snapshots  developed  for  Hurricane  Gilbert  Two  runs 
were  made.  The  first  generates  a  snapshot  wind  field  for  each  of  the  12 
Gilbert  cases  (including  realistic  forward  motion  and  steering  flow  parameters) 
and  interpolates  each  snapshot  to  a  polar  grid  using  a  test  history  table.  Inter¬ 
polations  are  also  made  from  pairs  of  adjacent  snapshots  (equal  time  weight). 
Winds  for  the  23  wind  fields  so  produced  on  the  polar  grid  are  then  azimuth- 
ally  averaged.  These  results  are  given  in  Appendix  C. 

A  second  test  run  was  made  which  modeled  Gilbert  during  its  passage 
across  the  Gulf  of  Mexico  between  1200  UT  15  September  through  0000  UT 
17  September,  1988.  Snapshots  for  this  run  consisted  of  Cases  7, 9, 11,  and 
12.  The  target  grid  for  this  run  was  a  grid  of  nominal  12  nm  spacing  covering 
the  Gulf  of  Mexico.  The  wind  fields  were  output  at  12-houriy  intervals. 
Additional  details  and  surface  wind  field  plots  are  given  in  Appendix  D. 
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Table  7 

Observed  Pressure  and  Azlmuthally  Averaged  Pseudo-Gradient  Wind  Maxima 
In  Hurricane  Gilbert1  and  Estimated  Generalized  Profile  Parameters 


Data/ 

Time 

(UTC) 

Input 

Output  | 

n 

g 

K- 

m/s 

Wi 

*P> 

mb 

mb 

B, 

■ 

1 

11/1624 

972 

1011 

47 

90 

40 

35 

47 

127 

23 

16 

2.52 

2.52 

2 

11/2002 

968 

1011 

58 

106 

47 

36 

58 

162 

32 

11 

2.52 

2.52 

3 

13/1744 

905 

1012 

16 

111 

65 

32 

16 

125 

94 

13 

1.68 

2.39 

4 

13/2333 

888 

1012 

a 

100 

69 

38 

13 

100 

110 

14 

1.59 

2.52 

S 

14/0649 

m 

1012 

13 

70 

62 

44 

13 

70 

101 

18 

1.41 

2.52 

6 

14/1126 

890 

1012 

13 

69 

62 

48 

13 

n 

102 

20 

1.41 

2.52 

7 

16/1204 

951 

1012 

a 

65 

25 

38 

22 

87 

45 

16 

0.56 

2.52 

8 

16/1629 

950 

1012 

32 

69 

28 

41 

32 

103 

35 

27 

0.94 

2.51 

8 

16/0016 

949 

1010 

El 

57 

27 

40 

81 

41 

El 

CO 

2.52 

10 

16/0539 

950 

1011 

55 

120 

41 

37 

55 

170 

52 

n 

1.12 

2.24 

11 

16/1850 

953 

1010 

50 

100 

42 

43 

50 

119 

47 

10 

1.33 

2.52 

12 

16/2131 

954 

1009 

a 

100 

39 

40 

40 

112 

45 

10 

1.19 

2.52 

|  '  Cases  correspond  to  Black  and  Wl 

doughty's  (1992)  and] 

fait  of  aircraft  da 

tain  Hun 

ricane  Gilbert 
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Table  8 

Comparison  of  Measured  Flight-Level  Wind  Maxima  and  Fitted 
Gradient  Wind  Maxima  for  Double  Exponential  Pressure  Profile1 

Inner  Ring 

Outer  Ring 

Measured 

Fitted 

Fitted  | 

wm 

E3HEI!3H 

9 

Vt  at 

JJee^ 

km 

m/s 

km 

m/a 

km 

m/s 

km 

m/s 

1 

47 

40 

47 

39.6 

90 

35 

110 

37.1 

35.7 

2 

S3 

47 

58 

46.6 

108 

36 

108* 

35.7* 

— 

3 

16 

65 

16 

66.7 

111 

32 

105 

33.7 

33.6 

4 

13 

69 

13 

70.3 

100 

38 

85 

37.2 

36.3 

5 

13 

62 

13 

63.5 

70 

44 

59 

46.2 

45.4 

6 

13 

62 

13 

63.7 

69 

48 

59 

48.4 

47.3 

7 

22 

25 

22 

26.3 

65 

38 

82 

39.2 

36.1 

rm 

32 

28 

30 

29.7 

69 

41 

98 

48.1 

39.2 

9 

22 

27 

21 

28.9 

57 

40 

78 

43.4 

38.1 

10 

55 

41 

55 

39.2 

120 

37 

120 

37.1 

— 

11 

ESI 

42 

Did 

100 

43 

93 

41.5 

41.3 

12 

40 

39 

ES 

32.9 

100 

40 

92 

38.4 

38.2 

1  Cases  correspond  to  Black  and  Willoughby's  (1992)  analysis  of  Hurricane  Gilbert 
*  Second  ring  maximum  not  resolved,  profile  gradually  decays  from  inner  maximum 
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Tangential  Wind  (m/s)  TangeniM  WM  (mtt)  TangenM  Wind  (mfi) 


Figure  5.  Comparison  of  azimuthaHy  averaged  reconnaissance  winds  and  fitted  gradient 
winds  in  12  cases  of  Hurricane  Gilbert  defined  by  Black  and  Willoughby  (1992) 
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4  Summary 


The  CE  tropical  cyclone  surface  wind  field  model  has  been  a  very  useful 
tool  in  ocean  response  modeling  for  more  than  a  decade.  The  model  continues 
to  be  used  regularly.  The  CE  recently  held  a  workshop  to  reassess  model 
assumptions,  particularly  in  light  of  modem  advances  in  computing  technology 
and  field  measurement  of  hurricane  structure.  Model  limitations  were  identi¬ 
fied  and  evaluated  in  terms  of  their  perceived  importance  to  ocean  response 
modeling  and  the  level  of  effort  required  to  develop  improved  solutions.  The 
limitations  are  summarized  in  this  report 

Two  aspects  of  the  CE  model  were  targeted  for  improvement.  This  report 
describes  the  improvements  developed  for  die  upgraded  model.  First  the 
standard  CE  model  represents  a  compromise  between  spatial  resolution  in  the 
central  region  of  very  high  gradients,  coverage  of  the  foil  ocean  area  affected 
by  the  tropical  cyclone,  and  computer  requirements.  Computing  resources  are 
much  more  available  now  than  at  the  time  the  model  was  developed  in  the  late 
1970’s.  The  model  was  upgraded  to  include  more  computationally  intensive 
options  which  give  improved  resolution  and  areal  coverage.  Up  to  seven 
nested  grids  are  now  available,  compared  to  only  five  nests  in  the  standard 
model.  In  a  typical  application,  this  upgrade  can  be  used  to  achieve  2-km 
resolution  around  foe  eye  (as  compared  to  5-km  resolution  often  used  in  the 
standard  model)  and  an  expanded  total  coverage  area. 

The  second  upgrade  to  the  standard  CE  model  allows  a  more  general  speci¬ 
fication  of  the  axisymmetric  pressure  profile.  This  upgrade  can  be  used  to 
create  wind  fields  with  maxima  at  two  different  radii  or  with  a  broad  maxi¬ 
mum  extending  over  a  range  of  radii.  It  also  provides  more  flexibility  in  fit¬ 
ting  the  shape  of  single  peaked  wind  profiles. 

The  upgraded  model  is  demonstrated  with  historical  hurricanes.  The  five- 
nest  and  seven-nest  models  are  applied  to  Hurricane  Camille.  The  fully 
upgraded  model,  with  seven  nests  and  general  pressure  specification,  is  applied 
to  Hurricane  Gilbert  This  hurricane  was  chosen  because  it  is  well-documented 
by  Black  and  Willoughby  (1992)  and  it  evolved  into  some  non-traditional 
storm  structures.  The  upgraded  model  was  more  effective  than  foe  standard 
CE  model  in  simulating  foe  storm. 
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Appendix  A 

Comparison  of  Five-Nest  and 
Seven-Nest  Models  for 
Hurricane  Camille 


Comparison  of  modeled  winds  at  20-m  height  and  at  30-min  intervals  in 
alternative  hindcasts  of  Hurricane  Camille  in  die  Gulf  of  Mexico  during 
August  1969.  The  run  labelled  “5-Nests”  used  the  standard  CE  model  with 
5 -km  spacing  on  the  inner  nest  The  run  labelled  “7-Nests"  used  the  upgraded 
CE  model  with  seven-nests  and  grid  spacing  of  2  km  on  the  inner  nest.  The 
first  (second)  line  at  each  time  step  gives:  maximum  scalar  (vector  magnitude) 
wind  speed  difference  found  on  the  grid  between  the  two  runs,  latitude  and 
longitude  of  die  grid  point,  and  wind  speed  and  direction  of  the  alternative 
solutions  at  that  grid  point. 
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A1 


Table  Al 

Comparison  of  Wind  Estimates  from  5-Nest  and 
Hurricane  Camille 


Orid  Point  Coord. 


7-Nest  Models, 


161330  2.2 

7.2 


161400  1.5 


161430  3.4 

7.3 


§4test  Modal 

74last  Modd 

sr- 

Olr. 

dag  as 

■t/a 

Dir. 

dag  as 

161200 

7.1 

24.1 

-85.7 

27.9 

151.7 

20.8 

162.6 

8.4 

24.1 

•85.7 

27.9 

151.7 

20.8 

162.6 

161230 

3.0 

24.1 

-86.0 

39.7 

336.6 

36.7 

342.5 

6.0 

24.1 

-85.7 

42.0 

172.2 

43.1 

180.1 

161300 

1.2 

24.1 

-85.7 

41.6 

181.0 

42.8 

183.7 

7.0 

24.3 

-86.0 

52.6 

47.1 

52.9 

54.6 

161500 


161530  2.4 

8.2 


161600  1.0 


161630  3.8 

6.2 


161700  1.5 


161730 


161800  2.0 
6.9 


161830  3.4 

7.0 


161900  2.1 


161930  2.8 

6.0 


162030 


(Shoot  1  oft) 
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Grid 


162100  1.2 
3.2 


162130  2.4 

6.8 


Point  Coord. 


54lnt  Modd 

7-Nest  Model 

n  . 

Mr. 

fj  . 

Mr. 

m/e 

dagae 

m/s 

dag  az 

44.1 

1S9.6 

45.3 

163.4 

44.1 

150.6 

45.3 

163.4 

13.3 

42.6 

132.4 

46.9 

170000 


170030 


170100  2.2 

7.2 


170130 

2.4 

25.6 

-87.4 

31.0 

5.0 

25.6 

-87.1 

45.7 

170200  0.9 

7.8 


170230 

7.3 

25.8 

-67.4 

19.8 

104.9 

iae 

124.3 

9.0 

25.8 

•87.4 

19.8 

104.9 

12.6 

124.3 

170300 

3.5 

25.8 

-87.4 

28.9 

221.5 

32.3 

226.8 

6.7 

26.0 

-87.4 

51.2 

81.4 

52.0 

88.8 

170330 

1.0 

26.0 

-87.6 

51.3 

7.4 

50.4 

14.3 

7.5 

26.0 

-87.4 

42.1 

121.2 

41.2 

131.4 

170400 

2.8 

26.0 

-87.6 

38.2 

333.4 

35.4 

336.4 

5.5 

26.0 

-87.4 

41.5 

175.3 

42.9 

182.5 

170430 

2.9 

26.2 

-87.6 

45.2 

39.4 

42.3 

48.6 

7.6 

26.2 

-87.6 

45.2 

39.4 

42.3 

48.6  1 

170500 

1.3 

26.2 

-87.4 

44.2 

161.7 

45.5 

164.7  1 

3.7 

26.4 

-87.6 

53.4 

53.5 

54.4 

57.3  | 

170530 

1.7 

26.2 

-87.6 

32.2 

250.2 

34.1 

251.5  1 

8.1 

26.4 

-87.6 

46.2 

78.9 

45.2 

89.0  1 
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QrM 


170800  6.0 

8.1 


170630  2.5 

5.2 


170700  4.0 

6.7 


170730  2.4 

3.4 


170600  3.2 

7.9 


170000  2.0 

6.3 


170030  0.9 


171000  5.0 

6.6 


171030  1.4 

6.4 


171100  2.0 

7.1 


171130  2.3 

5.3 


171200 


171230  3.7 

5.8 


171300  1.8 

7.1 


Point  Coord. 


MoM 

7-Maat  Modal 

,  . 

Dir. 

Ur. 

m/a 

dOBK 

m/a 

dag  ax 

22.5 

163.0 

15.7 

177.4 

22.5 

163.0 

15.7 

177.4 

37.7 

102.0 

40.2 

180.6 

37.7 

102.9 

40.2 

100.6 

42.8 

20.2 

38.8 

36.7 

42.8 

20.2 

38.8 

36.7 

6.5 

353.2 

4.1 

20.6 

6.5 

353.2 

4.1 

20.6 

30.1 

220.2 

33.2 

234.5 

51.4 

81.1 

51.1 

80.0 

40.3 

7.7 

47.8 

16.1 

41.8 

120.6 

40.9 

130.7 

35.2 

344.8 

32.3 

345.8 

43.0 

163.9 

45.6 

171.6 

8.8 

290.2 

8.0 

206.4 

49.4 

53.7 

48.5 

63.7 

342.6 

43.9 

186.0 

33.3 

171330 

1.7 

27.4 

-88.5 

40.1 

313.1 

38.4 

316.2 

7.5 

27.6 

-68.3 

50.7 

111.3 

51.0 

110.7 

171400 

1.4 

27.6 

-88.5 

28.0 

12.5 

26.6 

6.4 

5.0 

27.6 

-68.3 

48.2 

140.0 

40.5 

155.8 

171430 

2.3 

27.6 

-88.5 

30.2 

280.3 

32.6 

279.7 

8.4 

27.8 

-88.5 

48.7 

58.4 

47.3 

68.2 
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QrM  Point  Coord. 

6-Neat  Model 

7-Nest  Model  j 

Max 

Day/ 

Spaed 

Let 

Long. 

3pff(| 

Dir. 

Speed 

Oir. 

Hour 

DM.  an/s 

deg 

deg 

m/s 

deg  ax 

m/a 

deg  ax 

171 S00 

6.0 

27.8 

-88.5 

19.3 

114.2 

13.3 

125.6 

6.8 

27.8 

-88.5 

19.3 

114.2 

13.3 

125.6 

171530 

1.7 

27.8 

-88.7 

45.9 

338.5 

44.2 

342.8 

6.5 

28.0 

-88.5 

52.5 

61.4 

53.2 

88.4 

171600 

1.6 

28.0 

-88.5 

41.8 

126.4 

40.1 

136.8 

7.6 

28.0 

-88.5 

41.8 

126.4 

40.1 

136.8 

171630 

1.7 

28.0 

-88.7 

34.0 

338.7 

32.3 

338.4 

7.5 

28.0 

-88.5 

45.2 

169.1 

46.3 

178.4 

171700 

3.1 

28.3 

-88.7 

44.9 

45.2 

41.8 

54.2 

7.5 

28.3 

-88.7 

44.9 

45.2 

41.8 

54.2 

171730 

2.0 

28.3 

-88.7 

6.9 

324.1 

4.9 

340.6 

3.2 

28.5 

-88.7 

53.8 

59.8 

53.9 

63.2 

171800 

2.6 

28.3 

-88.7 

33.7 

240.9 

36.4 

246.3 

8.5 

28.5 

-88.7 

47.2 

81.6 

45.0 

91.9 

171830 

7.8 

28.5 

-88.7 

24.5 

167.5 

16.7 

175.6 

8.3 

28.5 

-88.7 

24.5 

167.5 

16.7 

175.6 

171900 

2.5 

28.5 

-88.7 

38.2 

214.3 

40.7 

221.1 

8.9 

28.7 

-88.7 

50.1 

97.5 

49.6 

107.7 

171930 

2.9 

28.7 

-89.0 

41.2 

1.6 

38.3 

6.5 

7.7 

28.7 

-88.7 

41.2 

154.8 

40.8 

165.5 

172000 

1.8 

28.7 

-88.7 

43.2 

192.8 

45.0 

197.7 

7.2 

28.9 

-89.0 

50.0 

38.1 

49.0 

46.4 

172030 

1.5 

28.9 

-88.7 

47.9 

150.2 

49.4 

158.3 

7.1 

28.9 

-88.7 

47.9 

150.2 

49.4 

158.3 

172100 

2.0 

28.9 

-89.0 

31.1 

283.6 

33.2 

283.5 

8.2 

29.1 

•89.0 

47.6 

55.6 

46.1 

65.4 

172130 

3.1 

29.1 

-89.0 

15.6 

79.7 

12.5 

86.5 

3.5 

29.1 

-89.0 

15.6 

79.7 

12.5 

86.5 

172200 

3.6 

29.1 

-89.0 

27.7 

250.0 

31.3 

255.4 

8.5 

29.3 

-89.0 

50.1 

72.7 

49.6 

82.6 

172230 

7.2 

29.3 

-89.0 

23.8 

135.6 

16.6 

147.2 

8.2 

29.3 

-89.0 

23.8 

135.6 

16.6 

147.2 

172300 

1.6 

29.3 

•89.0 

37.8 

205.3 

39.4 

214.6 

7.1 

29.5 

-89.0 

51.5 

93.0 

52.0 

100.9 

172330 

2.2 

29.5 

-89.2 

44.9 

7.8 

42.7 

15.1 

6.2 

29.5 

-89.0 

39.1 

143.3 

37.6 

152.3 
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omnBMHi 

mm 

arid  Point  Coord. 

S  Neal 

Modal 

7-Neat  Model  || 

Day/ 

fc'isii 

Lat 

Long. 

Dir. 

Sptod 

Dir. 

Hour 

*9 

*9 

E2fll 

ml* 

deg  sz 

180000 

1.7 

29.5 

-89.0 

43.4 

188.7 

45.1 

195.4 

5.5 

20.5 

-89.0 

43.4 

188.7 

45.1 

195.4 

1.8 

20.7 

•69.2 

27.4 

30.7 

25.6 

26.6 

6.0 

28.7 

-89.0 

49.7 

142.6 

50.7 

149.4 

180100 

3.2 

28.7 

-89.2 

27.0 

274.3 

30.2 

276.4 

8.4 

28.0 

-80.2 

50.6 

62.6 

49.5 

72.1 

180130 

6.0 

28.0 

-88.2 

20.1 

126.8 

13.2 

141.7 

8.0 

20.0 

-89.2 

20.1 

126.8 

13.2 

141.7 

180200 

1.4 

29.0 

-89.2 

36.5 

197.3 

37.9 

206.1 

7.1 

28.8 

-89.2 

36.5 

197.3 

37.9 

208.1 

180230 

2.0 

30.1 

•80.4 

43.9 

16.8 

41.0 

24.2 

0.2 

30.1 

•89.2 

46.0 

132.5 

45.5 

143.9 

1.4 

30.1 

-80.2 

45.9 

176.3 

47.3 

181.7 

BUM 

6.0 

30.3 

-80.4 

52.0 

47.1 

52.1 

53.6 

180330 

2.5 

30.3 

-80.4 

23.1 

71.2 

20.6 

66.5 

3.3 

30.3 

-89.2 

49.9 

130.7 

50.5 

143.4 

180400 

1.0 

30.3 

-89.4 

24.1 

245.7 

26.0 

252.7 

3.6 

30.3 

89.4 

24.1 

245.7 

26.0 

252.7 

180430 

0.0 

28.3 

82.2 

6.4 

293.0 

7.5 

200.6 

1.4 

30.3 

89.4 

38.9 

234.0 

30.7 

235.6 

180500 

0.0 

28.3 

82.2 

8.3 

201.2 

7.5 

296.5 

1.4 

20.5 

84.0 

4.4 

353.8 

4.3 

11.9 

E 

0.0 

28.5 

82.2 

8.8 

202.6 

7.9 

296.1 

mm 

1.5 

30.3 

88.4 

39.5 

236.7 

38.8 

234.8 

180600 

0.8 

28.5 

82.4 

7.8 

293.4 

7.0 

299.2 

1.4 

20.5 

83.3 

6.5 

334.2 

6.0 

345.8 

180630 

0.8 

28.7 

82.4 

7.5 

203.6 

6.6 

200.3 

1.3 

29.7 

82.9 

8.0 

328.1 

7.2 

336.0 

180700 

0.0 

29.5 

82.6 

7.6 

315.2 

6.7 

322.1 

1.3 

29.5 

82.6 

7.6 

315.2 

6.7 

322.1 

180730 

0.8 

29.5 

82.6 

6.7 

310.7 

5.8 

317.2 

1.1 

29.5 

82.6 

6.7 

310.7 

5.8 

317.2 

0.8 

29.5 

82.4 

6.8 

301.6 

5.9 

306.7 

1.1 

29.5 

82.4 

6.8 

301.6 

5.9 

306.7 

180830 

0.0 

29.7 

82.2 

7.2 

298.6 

6.3 

302.9 

1.0 

29.7 

82.0 

8.2 

294.2 

7.3 

298.0 
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Appendix  B 

Documentation  of  CE  Model 
Upgrades 


This  appendix  provides  brief  documentation  of  new  and  modified  programs 
in  the  upgraded  CE  tropical  cyclone  surface  wind  field  model.  The  material  is 
a  supplement  to  the  primary  documentation  of  Caidone  et  al.  (1992).  Five 
FORTRAN  programs  are  discussed.  The  programs  HIST_ADC.7NE, 
SNAP_ADC.7NE,  and  SNAP_HOL.7NE  are  modified  versions  of  die  previous 
HIST  and  SNAP  programs.  Programs  1EYEWALL.HOL  and 
2EYEWALL.HOL  are  new.  They  are  helpful  in  implementing  the  new  gener¬ 
alized  surface  pressure  specification. 


Program  HIST_ADC.7NE 

HIST_ADC.7NE  is  a  slight  modification  of  HIST_ADC.F.  All  input  files 
except  LSNAP,  and  all  output  files,  are  unchanged  from  HIST_ADC.F  to 
HIST_ADC.7NE.  The  changes  to  file  LSNAP  are  as  follows: 

-  pressure  arrays,  formerly  dimensioned  (21,21,5),  are  now  dimensioned 
(21,21,7); 

-  wind  arrays,  formerly  dimensioned  (21,21,10),  are  now  dimensioned 
(21,21,14); 

-  in  the  variables  at  die  aid  of  records  in  LSNAP,  variable  DX  is  fol¬ 
lowed  by  the  new  variable  INSIDE.  INSIDE  is  an  integer  indexing  the 
innermost  live  nest  of  the  7  nests  supported,  so  that  the  effective  grid 
spacing  is  DX*2**(INSIDE-1) . 


Program  SNAP_ADC.7NE 

SNAP_ADC.'TNE  is  derived  from  SNAP_ADC.F.  The  modifications  allow 
the  user  to  run  up  to  7  grid  nests  rather  than  the  previous  mandatory  5  nests. 
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B1 


Arrays  of  pressure,  pressure  gradient,  wind,  formerly  dimensioned  (21,21,5), 
are  now  dimensioned  (21,21,7).  Two  changes  are  made  to  namelist  /NAMES/ 
as  follows: 

-  a  variable  name  NOMEN  (CHARACTERS)  for  storm  identification 
has  been  included  ; 

-  integer  variable  INSIDE  has  been  added.  INSIDE  indexes  the  finest 
live  nest  of  the  7  nests  provided.  Thus  the  number  of  live  nests  is 
(8- INSIDE),  and  the  grid  spacing  of  foe  finest  live  nest  is 
DX*2**0NSIDE-1).  Default  values  are  DX  =  2.  and  INSIDE  =  1, 
yielding  a  2  km  grid  spacing  and  an  execution  time  roughly  4  times  as 
long  as  the  existing  5-nest  model.  The  combination  DX  =  6.25, 
INSIDE  =  3,  (lines  304,  305)  reproduces  the  25  km  spacing  often  used 
by  the  CE  for  global  studies.  In  the  great  majority  of  applications,  the 
useful  values  of  INSIDE  are  1,  2,  and  3.  INSIDE  =  4  may  be  tried  for 
running  a  quick  preliminary  study  on  a  coarse  grid. 


Program  SNAPHOL7NE 

SNAP_HOL.7NE  is  an  extensive  modification  of  SNAP_ADC.7NE  to 
include  the  generalized  pressure  profile  as  well  as  foe  capability  for  modeling 
up  to  7  nests.  The  variable  ST12  and  tire  quadrantal  variation  of  PFAR  and 
RADIUS  have  been  excised.  In  OWI’s  experience  with  foe  hurricane  model, 
they  have  been  used  only  once:  for  hurricane  Efoise,  September  13,  1975. 
Variable  ITRACK  has  been  excised:  its  use  pertained  to  a  1969  study  in  which 
direction  was  specified  in  points.  Namelist  /NAMES/  is  changed  materially. 
Each  variable  and  array  in  /NAME3/  is  documented  below.  The  method  of 
computation  of  pressure  and  pressure  gradient  is  discussed  in  a  later  part  on 
Mathematical  Method. 


Revisions  at  beginning  of  program 

REAL  RADIUS(2),DPRESS(2),HOLL(2) 

CHARACTERS  NOMEN 

EQUIVALENCE  (RADI RADIUS),  (RAD2RADIUS(2)),  (Bl.HOLL), 
$  (B2ROLL(2)),  (DPI, EXPRESS),  (DP2,DPRESS(2)) 

NAMELIST  /NAME3/  SGW,  AN1,  NOMEN, 

$  EYELAT,  EYLONG,  DIREC,  SPEED,  EYPRES,  RADIUS,  RADI, 

$  RAD2,  PFAR,  NM,  DX,  INSIDE,  HOLL,  Bl,  B2,  DPRESS,  DPI, 

$  DP2 


B2 
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Definition  of  variables  In  namelist  NAME3 

SGW  Magnitude  of  surface  geostrophic  wind,  m/sec 
AN1  Angle  between  SGW  and  east,  counterclockwise  from  east 
NOMEN  Designator  for  tropical  storm,  e.g.  two  digits  and  one  letter 
EYELAT  Latitude  of  eye  of  storm  at  snap  tune,  north  positive 
EYLONG  Longitude  of  eye  of  storm  at  snap  time,  east  positive  (EYLONG  is 
included  for  archival  purposes  but  not  presently  used  in 
computation) 

DIREC  Direction  of  forward  motion  of  storm,  clockwise  from  north 

SPEED  Speed  of  forward  motion,  in  kt  (but  redefinable  according  to 

die  switch  variable  UNITS) 

EYPRES  Pressure  at  eye  of  storm,  in  mb 

RADIUS  Scale  radius  of  the  two  components  of  exponential  pressure  profile 
RAD1JLAD2  Alternate  names  for  specifying  RADIUS;  convenient  when  only 
one  exponential  is  modeled 
PFAR  Ambient  pressure  exterior  to  storm,  in  mb 

NM  Number  of  computational  cycles  in  nest  1;  NM  should  be  a 

multiple  of  64  (default:  NM  =  3200) 

DX  Grid  spacing  in  nest  1,  in  km  (default:  DX  =  2.) 

INSIDE  Index  of  finest  nest  actually  used  for  computations;  the  finest 

active  grid  spacing  is  DX*2**(INSIDE-1)  (default:  INSIDE  =  1) 
HOLL  Power  to  which  radius  is  raised  in  die  modified  Holland’s  (1980) 
pressure  profile  model.  When  HOLL(2)  =  0.,  only  one  exponen¬ 
tial  is  used;  RADIUS(2)  and  DFRESS(2)  are  then  ignored. 

B1.B2  Alternate  names  for  specifying  HOLL;  convenient  when  only  one 
exponential  is  modeled.  Default  values  are  B1  =  1.,  B2  =  0. 

These  defaults  are  reinstated  for  every  snapshot.  Use  of  the 
defaults  reverts  to  a  standard  exponential  pressure  profile,  as  used 
in  SNAP_ADC.7NE. 

DPRESS  DPRESS©  is  the  pressure  difference  associated  with  RADIUS(I) 
and  HOLL(I)  in  OWI’s  double -eye  wall  extension  of  Holland’s 
(1980)  modified  exponential  profile. 

DPI  J5P2  Alternate  names  for  specifying  DPRESS.  As  explained  in  the  part 
on  Mathematical  Method,  it  is  never  advantageous  to  include  DP2 
in  an  input  list;  it  is  in  die  NAMELIST  in  order  to  force  its 
appearance  in  die  output  file. 

SPECIAL  WARNING:  Do  not  input  both  members  of  an  equivalence.  Input 
either  RADIUS  or  RAD  and  RAD2;  either  HOLL  or  B1  and  B2;  either 
DPRESS  or  DPI  and  DP2.  The  program  imposes  consistency  checks  on 
RADIUS,  HOLL,  and  DPRESS;  it  does  not  check  EYPRES  and  PFAR,  so  that 
it  remains  the  user’s  duty  to  verify  that  PFAR  >  EYPRES. 
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Mathematical  method 


B4 


The  program  performs  the  foilowing  consistency  checks: 

1.  If  B1  <  0  or  B2  <  0,  stop. 

2.  If  B2  =  0,  run  Holland’s  modified  exponential  model  (ignore  RAD2 
and  DP2). 

2.1  If  RADI  £  0,  stop. 

2.2  If  DPI  not  specified,  compute  it  as  PFAR-EYPRES. 

2.3  If  DPI  specified,  but  inconsistent  with  PFAR-EYPRES,  stop. 

2.4  Set  N  =  1. 

3.  If  B2  >  0,  run  OWI’s  extension  of  Holland's  pressure  profile. 

3.1  If  RADI  £  0  or  RAD2  £  0,  stop. 

3.2  If  DP2  not  specified,  compute  it  as  PFAR-EYPRES- DPI . 

3.3  If  DP2  specified,  but  inconsistent  with  PFAR-EYPRES-DP1 ,  stop. 

3.4  Set  N  =  2. 


Program  1  EYEWALL.HOL 

Program  1  EYEWALL.HOL  attempts  to  fit  snapshot  parameters  to  a  guessed 
wind  profile.  It  considers  only  die  case  of  one  exponential  (B2  =  0.0). 
1EYEWALL.HOL  requires  the  following  input  arguments  On  namelist  /INN/): 

BLAT  Absolute  value  of  latitude  of  eye,  degrees  &  decimals;  used  in 

computation  of  die  Coriolis  parameter.  Also,  die  value  BLAT=99 
is  used  as  a  flag  to  stop  computation. 

EYFRES  Pressure  at  eye,  in  mb 

PFAR  Pressure  at  large  (theoretically  infinite)  distance  from  eye,  in  mb 
RW1  Radius  at  which  a  wind  speed  is  guessed,  in  km 
SP1  Wind  speed  corresponding  to  RW1,  in  m/sec 
V150  Wind  speed  at  radius  of  150  km,  in  m/sec 

The  following  outputs  are  primed: 

1.  In  namelist  /INN/: 

SP11,  SP12  a  SP1  minus  &  plus  1  m/sec 

SP31,  SP32  a  Y150  minus  &  plus  1  m/sec 

DP2  s  pressure  difference  (far  field  minus  eye),  in  pascal 

2.  In  namelist  /VORTEX/: 

COR  a  Coriolis  parameter,  2*omega*sin(BLAT) 
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FR22  =  the  quantity  0.5*COR*RW1  (used  in  the  computation  of 
gradient  wind) 

FR23  =  0.5*COR*r,  where  r  is  150000  m  or  150  km 
re  YE  =  pressure  at  eye,  in  pascal 

3.  Below  namelist  /VORTEX/,  six  parameters  are  printed,  defined  from  left  to 
right  as: 

3.1  DP2  (see  above) 

3.2  Fitted  value  of  scale  radius,  in  m 

3.3  Fitted  value  of  Holland’s  exponent 

3.4  Fitted  value  of  gradient  wind  at  radius  RW1,  in  m/sec 
(in  a  good  fit,  this  will  be  neatly  equal  to  SP1) 

3.5  Fitted  value  of  gradient  wind  at  radius  150  km,  in  m/sec 
(in  a  good  fit,  this  will  be  nearly  equal  to  V150) 

3.6  Goodness  of  fit  measure:  absolute  value  of  gradient  wind  minus 
in  put  wind  at  radius  RW1,  plus  the  same  at  radius  150  km.  A  value 
less  than  3.0  implies  a  tolerably  well-fitting  solution. 

4.  Table  with  6  columns  and  150  lines: 

4.1  Radius,  in  km 

4.2  Pressure,  in  mb 

4.3  First  component  of  pressure  gradient,  in  pascal/m 

4.4  Second  component  (this  is  zero,  since  only  one  exponential  was 
fitted) 

4.5  Pressure  gradient  (here  equal  to  output  #4.3) 

4.6  Gradient  wind,  m/sec 

5.  Below  the  table  are  numbers  that,  if  the  fit  is  satisfactory,  the  user  can 
insert  into  namelist  /NAME3/  of  SNAP_HOL.7NE: 

EYELAT  s  echo  of  the  input  BLAT 
EYPRES  =  echo  of  input 
PFAR  =  echo  of  input 
RADI  =  scale  radius,  in  nm 

HOLL  =  two  values  of  Holland’s  exponent;  the  second  value  is  zero, 
because  oily  one  exponential  was  fitted 


Program  2EYEWALL.HOL 

Program  2EYEWALL.HOL  attempts  to  fit  snapshot  parameters  to  a  guessed 
wind  profile.  It  considers  the  case  of  two  exponentials  (B2  >  0.0). 
2EYEWALL.HOL  requires  the  following  inputs  in  namelist  /INN/: 

BLAT  Same  usage  as  in  1EYEWALL.HOL 
EYPRES  Same  usage  as  in  1EYEWALL.HOL 
PFAR  Same  usage  as  in  1 E YE W ALL.HOL 
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RSI  Scale  radius  of  inner  ring,  in  km 

(numerical  experiments  with  this  scheme  have  shown  that  the  inner 
scale  radius  can  safely  be  taken  equal  to  die  inner  radius  to  local 
maximum  wind) 

RW2  Radius  to  maximum  wind  of  outer  ring,  in  km 

DRING  An  integer  switch  variable  indexing  die  shape  of  die  wind  profile 

DRING  =  1:  the  maximum  wind  is  greater  in  the  inner  ring 
DRING  *  2:  the  maximum  wind  is  greater  in  the  outer  ring 

SP1  Desired  wind  speed  at  radius  RSI,  in  m/sec 

SP2  Desired  wind  speed  at  radius  RW2,  in  m/sec 

The  following  outputs  are  printed: 

1.  In  namelist /INN/: 

SP11,  SP12  =  SP1  minus  and  plus  1  m/sec 
SP21,  SP22  =  SP2  minus  and  plus  1  m/sec 

2.  In  namelist  /VORTEX/: 

BLAT  =  echo  of  input 
COR  =  same  usage  as  in  1EYEWALL.HOL 
RADI  =  RSI,  in  m 
RAD2  =  RW2,  in  m 

RAD3  =  0.5*(RAD1+RAD2);  a  local  minimum  of  wind  speed,  if 
found,  will  be  near  RAD3 
RAD4  =  130000  m  (=  130  km) 

FR21  =  die  quantity  0.5*COR*RAD1;  used  in  the  computation  of 
gradient  wind 

FR22  =  the  quantity  0.5*COR*RAD2 
FR23  =  the  quantity  0.5*COR*RAD3 
PEYE  =  same  usage  as  in  1EYEWALL.HOL 
DP  =  pressure  difference  (far  field  minus  eye),  in  pascal 

3.  Below  namelist  /VORTEX/,  nine  parameters  are  primed,  defined  from  left 
to  right  as: 

3.1  Fitted  value  of  DPI  (partial  pressure  difference  for  first 
exponential),  in  pascal 

3.2  Fitted  value  of  DP2  (partial  pressure  difference  for  second 
exponential),  in  pascal 

3.3  Scale  radius  of  second  exponential,  in  m  (die  scale  radius  of  first 
exponential  has  been  fixed  at  RADI) 

3.4  Exponent  for  first  exponential  (Holland  1980) 

3.3  Exponent  for  second  exponential  (OWI  extension  of  Holland  (1980)) 

3.6  Fitted  wind  speed  at  radius  RADI,  m/sec 

3.7  Fitted  wind  speed  at  radius  RAD2,  m/sec 
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3.8  Fined  wind  speed  at  radius  RAD3.  m/scc 

3.9  (printed  below  3.1):  fitted  wind  speed  at  radius  130  km 

4.  Second  printing  of  namelist  /INN/,  if  given: 

Hie  fit  in  the  above  two-line  summary  was  unsatisfactory  in  that  the  wind 
at  150  km  was  greater  than  the  wind  at  RAD2;  a  second  fit  will  be 
attempted,  this  time  minimizing  the  wind  speed  at  130  km. 

SP11,  SP12  =  SP1  minus  and  plus  2  m/sec 
SP21,  SP22  =  SP2  minus  and  plus  2  m/sec 

5.  Below  second  printing  of  /INN/:  the  same  nine  parameters  as  in  #3  above, 
for  the  second  attempted  fit 

6.  Table  with  6  columns  and  130  lines: 

6.1  Radius,  in  km 

6.2  Pressure,  in  mb 

6.3  First  component  of  pressure  gradient,  in  pascal/m 

6.4  Second  component  of  pressure  gradient,  in  pascal/m 

6.5  Pressure  gradient  (sum  of  the  two  components),  in  pascalAn 

6.6  Gradient  wind,  in  m/sec 

7.  Below  the  table  are  numbers  that,  if  the  fit  is  satisfactory,  the  user  can 
insert  into  namelist  /NAME3/  of  SNAPJHOL.7NE: 


EYELAT 

EYPRES 

PFAR 

RADIUS 

HOLL 

DPI 


echo  of  die  input  BLAT 
echo  of  input 
echo  of  input 

two  values  of  scale  radius,  in  nm 

two  values  of  Holland’s  exponent 

pressure  difference  for  first  component,  in  mb 

(program  SNAP.HOL.7NE  computes  DP2  by  subtraction) 
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Appendix  C 
Sample  Application  of 
Upgraded  CE  Model  to  Simu¬ 
lation  of  12  Snapshots  of 
Hurricane  Gilbert 


This  appendix  provides  input  file  information  used  by  OWI  in  test  runs 
with  the  upgraded  CE  model,  including  seven  nests  and  the  generalized  pres¬ 
sure  specification.  Snapshots  were  generated  at  6-hr  intervals  for  Hurricane 
Gilbert,  which  occurred  during  September  1988.  The  first  snapshot  represents 
1200  UTC  (Universal  Tune  Coordinate,  formerly  known  as  Greenwich  Mean 
Tune)  IS  September  1988.  In  all,  12  snapshots  were  generated.  Inputs  for  the 
programs  SNAP.HOL.7NE  and  HIST.ADC.7NE  arc  included  as  implemented 
on  die  OWI  VAX  computer.  The  listed  values  of  parameter  AN1  follow  a 
meteorological  convention  (deg  azimuth  coming  from )  rather  than  the  conven¬ 
tion  used  by  die  Cardone  et  al.  (1992). 

Listings  of  die  full  field  of  surface  (19-m  elevation)  wind  speed  and  direc¬ 
tion  were  generated  on  a  polar  output  grid.  They  are  included  here  at  6-hr 
intervals  (snapshot  times).  Wind  speed  is  in  m/sec.  Wind  direction  is  in  deg 
azimuth  coming  from.  Printed  output  of  the  azimuthally  averaged,  surface 
wind  speed  and  inflow  angle  is  also  given  in  this  appendix.  In  addition  to  the 
12  snapshot  wind  fields,  this  output  includes  a  wind  field  interpolated  halfway 
betwee  each  pair  of  snapshots. 


Aaoandx  C  Smote  Aooficaion  12  Snmhote  of  Hurricano  Gabon 


Cl 


S  ASSIGN 

GRID. 312  P0R012 

S  ASSIGN 

05GILSERTZ1  F3R005 

S  ASSIGN 

GILB.WINOZ1  FOR052 

A  RUN  M0LL3 

SNAME1 

AZM 

8309, 

AON 

15XZ00, 

AMIN 

180, 

OX 

2.000030 

9 

ASTRES 

0, 

NSTRES 

17227, 

AW  I  NO 

X9, 

NWINO 

3X2, 

HH 

500.0030 

t 

INSIDE 

1, 

ATIHE 

1 

SENO 

SNAME2 

EYELAT 

16.00030 

OIREC 

290.0030 

*  Snapshot  1 

SPEED 

11.03030 

EYPRES 

972.0030 

PFAR 

1011.030 

RA01 

25.38030 

RA02 

68.73030l 

RADIUS 

25.38030 

,  68.73000  t 

OP1 

22.85030. 

OP2 

16.15030 

OPRES 

22.85030 

,  16.15000  , 

B1 

2.520030 

82 

2.520030 

9 

MOLL 

sew 

7.000030 

AN1 

110.0030 

ST12 

o.ooooosos+oo 

SENO 

SNAME2 

EYELAT 

16.00030 

OIREC 

290.0030 

*  Snapshot  2 

SPEEO 

11.00030 

EYPRES 

968.0030 

PFAR 

1011.030 

RADI 

31.32030 

RA02 

87.37030 

RAOIUS 

31.32030 

,  87.37000  , 

OP1 

31.77030 

0P2 

11.23030 

OPRES 

31.77030 

,  11.23000  , 

ai 

2.520030 

82 

2.520030 

HOLL 

2*2. 520030 

SGW 

7.000030 

AN1 

110.0030 

STI2 

0.0000030c*90 

SENO 

Figure  Cl.  Program  inputs,  Hurricane  GHbeit  (Sheet  1  erf  6) 
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C4 


Appendix  C  Sample  Appicaton,  12  Snapshots  of  Hurricane  Qflbert 


SNAHE2 

SYELA7 

x 

23.000)0 

OIREC 

SPEED 

x 

m 

290.00)0 

11.000)0 

*  Snapshot  9 

EYPRES 

9 

949.00)0 

PFAR 

9 

1010.0)0 

RA01 

9 

11.830)0 

RA02 

9 

A3. 530)0 

RADIUS 

9 

11.880)0 

,  43.53000  t 

DPI 

9 

40.670)0 

OP2 

X 

20.330)0 

DPRES 

X 

40.670)0 

20.33000  . 

81 

X 

0.75000)0 

62 

X 

2.5200)0 

HOLL 

X 

0.75000)0 

,  2.  520000  t 

SGW 

X 

7.0000)0 

AN1 

X 

110.00)0 

ST12 

SEND 

SNAHE2 

X 

0 .0009030c*  00 

* 

EYELA7 

X 

23.000)0 

OIREC 

X 

290.00)0 

•  Snapshot  10 

SPEED 

X 

11.000)0 

EYPRES 

X 

950.00)0 

PFAR 

X 

1011.0)0 

RADI 

X 

29.700)0 

RA02 

X 

91.630)0 

RADIUS 

X 

29.700)0 

f  91.63000  . 

DPI 

X 

51.840)0 

OP2 

X 

9.1600)0 

OPRES 

X 

51.840)0 

,  9.160000  , 

ai 

X 

1.1200)0 

B2 

X 

2.2400)0 

HOLL 

X 

1.1200)0 

l  2.240000  , 

SGW 

X 

7.0000)0 

AN1 

X 

110.00)0: 

ST12 

SEND 

SNAME2 

X 

0.000 00)0 E* 00 

EYELAT 

X 

24.000)0 

OIREC 

X 

290.00)0- 

,  Snapshot  1 1 

SPEED 

m 

11.00090 

EYPRES 

9 

953.00)0 

PFAR 

X 

1010.0)0 

RADI 

X 

27.000)0- 

RAD2 

X 

64.210)0 

RADIUS 

X 

27.000)0 

\  64.21000  » 

DPI 

X 

46.620)0 

DP2 

X 

10.380)0 

OPRES 

X 

46.62090 

10.38000  » 

BI 

X 

1.3300)0 

B2 

X 

2.5200)0 

HOLL 

X 

1.330090 

\  2.520000  * 

SGW 

X 

7.000090 

AN1 

X 

110.00)0 

S712 

SEND 

X 

0.0000030c*  00 

Figured.  (Sheet 4 of 6) 


12  Snapshots  of  Hurricane  Gifcert 


C5 


SNAME2 

EYELAT 

B 

24.00030  * 

OIREC 

B 

290.00301  » 

Snapshot  12 

SPEEO 

* 

11.00030  i 

EYPRES 

9 

954.0030  • 

P  FAR 

m 

1009.030  * 

RADI 

9 

21.60030  . 

RA02 

* 

60.61030  > 

RAOIUS 

at 

21.60030  * 

60.61000  > 

0P1 

M 

44.08030  • 

0P2 

9 

10.02030  , 

OPRE5 

* 

44.98030  , 

10.02000  * 

SI 

9 

1.190030  , 

a2 

9 

2.520030  , 

HOLL 

m 

1.190030  , 

2.  520000  . 

S  GW 

m 

7.000030  , 

AN1 

s 

110.0030  t 

ST12 

s 

G. 0000030c* 00 

SEND 

SNAME2 

EYELAT 

9 

999.0030  * 

OIREC 

9 

290.0030  i 

SPEEO 

9 

11.00030  » 

cYPRES 

9 

954.0030  , 

PEAR 

9 

1009.030  * 

RADI 

9 

0 • 000 0030 E* 00* 

RA02 

9 

0.00000306*00* 

RADIUS 

9 

2*0.00003006*00 

* 

OPl 

9 

0.00000306*00* 

DP2 

3 

0.00000106*00* 

OPRES 

M 

2*0.00003006*00 

9 

B1 

9 

1.000030i  » 

82 

9 

0.0000030E*  00* 

HOLL 

m 

1.000030  * 

0.00000006*00* 

sew 

• 

7.000010  * 

AMI 

9 

110.0030!  > 

ST12 

S 

0.00000306*00 

SEND 

0  20 

0  0  3  1 

0 

2  20 

0  0  3  2 

0 

A  20 

0  0  3  t  3 

0 

6  20 

0  0  3  •  4 

0 

8  20 

0  0  3  5 

0 

10  20 

0  0  16 

0 

12  20 

0  0  3  r  7 

0 

14  20 

0  0  3  8 

0 

16  20 

0  0  3  9 

0 

18  20 

0  0  3  10 

0 

20  20 

0  0  3  1  1 

0  . 

22  20 

0  0  3  1  2 

0 

999  0 

0  0  3  0 

0 

Figure  Cl .  (Sheet  5  of  6) 


Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gifcert 


HOUR  LAT 


0 

20 

0 

1 

20 

0 

2 

20 

0 

3 

20 

0 

4 

20 

0 

5 

20 

0 

6 

20 

0 

7 

20 

0 

8 

20 

0 

9 

20 

0 

10 

20 

0 

11 

20 

0 

12 

20 

0 

13 

20 

0 

14 

20 

0 

15 

20 

0 

16 

20 

0 

17 

20 

0 

18 

20 

0 

19 

20 

0 

20 

20 

0 

21 

20 

0 

22 

20 

0 

SWHAT 

KSTEP2 

S 

SEND 

workk: 

WORK*: 

WORKK: 

workk: 

workk: 

workk: 

workk: 

WORKK: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

workk: 

END  OF  STRESS 


STORM 

HISTORY 

L3NS 

SNAPS 

0 

0 

1  0 

0 

0 

1  2 

0 

0 

2  0 

0 

0 

2  3 

0 

0 

3  0 

0 

0 

3  4 

0 

0 

4  0 

0 

b 

4  5 

0 

0 

5  0 

0 

0 

5  6 

0 

0 

6  0 

0 

0 

6  7 

0 

0 

7  0 

0 

0 

7  8 

0 

0 

8  0 

0 

0 

8  9 

0 

0 

9  G 

0 

0 

9  10 

0 

0 

10  0 

0 

0 

10  11 

0 

G 

11  o 

0 

0 

11  12 

0 

0 

12  0 

22 

3 

880  9 

1 

880  9 

l 

880  9 

3 

880  9 

V 

880  9 

5 

88  0  9 

5 

880  9 

T 

880  9 

3 

880  9 

9 

860  9 

13 

880  9 

11 

880  9 

12 

880  9 

13 

880  9 

U 

880  9 

15 

880  9 

15 

880  9 

17 

880  9 

13 

880  9 

19 

860  9 

23 

880  9 

21 

880  9 

22 

880  9 

ZERO  HOUR  IS 
INTERP  ROT 


0.0000  0 

0.5000  0 

0.0000  0 
0.5000  0 

0.0000  9 

0.5000  9 

0.0000  0 

0.5000  9 

0.0000  9 

0.5000  9 

0.0000  0 
0.5000  0 

0.0000  9 

C. 5000  0 

0.0000  9 

0.5000  0 

0.0000  0 
0.5000  0 

0.0000  9 

0.5000  0 

0.0000  9 

C. 5000  9 

0.0000  9 


151200 

151509 

151809 

152100 

160000 

160309 

160600 

160900 

161209 

161509 

161800 

162109 

170009 

170300 

170609 

170900 

171200 

171500 

171800 

172100 

180000 

180309 

130609 


8809  151200 


Figured.  (Sheet 6 of 6} 


Appendx  C  Sample  Application,  12  Snapshots  of  Hurricane  Gifcert 


C7 


mi  >.««!  1*9450 


C8 


Appendx  C  Sample  Application,  12  Snapahota  of  Hurricane  Gtfbert 


Figure  C2.  Wind  speed  and  direction  fields,  Snapshot  1,  Hurricane  Gitoert  (Continued) 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gifeert 


Figure  C3.  Wind  speed  and  direction  fields.  Snapshot  2,  Hurricane  Glbert  (Continued) 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gifcert 


Figure  C4.  Wind  speed  and  direction  fields,  Snapshot  3,  Hurricane  Gibert  (Continued) 
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Appendx  C  Sample  Application,  12  Snapehoti  of  Hurricane  OHbert 


Figure  C6.  Wind  speed  and  direction  fields,  Snapshot  5,  Hurricane  GHbeit  (Continued) 
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Appendx  C  Simple  Application,  12  Snapshots  of  Hurricane  GMbert 


Figure  C7.  (Concluded) 
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Appendix  C  Sample  Application,  12  Snapshots  ot  Hurricane  Gitoert 


Figure  C8.  Wind  speed  and  direction  fields,  Snapshot  7,  Hurricane  Gitoert  (Continued) 
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Appendix  C  Sample  Application.  12  Snapshots  of  Hurricane  GHbert 
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Appendix  C  Sample  Application.  12  Snapshots  of  Hurricane  Gilbert 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gifeert 


Figure  Cl  1.  Wind  speed  and  direction  fields,  Snapshot  10,  Hurricane  Gibert  (Continued) 


AppwKfxC  Sampte  Appicalion,  12  SnapthoU  of  Hurricane  Ottart 


C27 


Figure  C11.  (Conducted) 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  GHbeit 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gilbert 


Figure  C13.  Wind  speed  and  direction  fields,  Snapshot  12,  Hurricane  Gibed  (Continued) 
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Figure  C13.  (Concluded) 
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Figure  C14.  Azimuthalty  averaged  speed  and  inflow  angle,  Hurricane  Gilbert  (Sheet  1  of  12) 
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Appendix  C  Sample  Application.  12  Snapshots  of  Hurricane  Gilbert 
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Figure  Cl  4.  (Sheet  3  of  12) 
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13. 5521 

29.3917 

180. 

333.360 

3.33  52 

7.4333 

17.3851 

24b. 

444.480 

7.  32  89 

4.0952 

12.1178 

360. 

666.720 

7.6295 

1. 4431 

11.4411 

480. 

888.960 

7.57  39 

0.6301 

11.7333 

600. 

1111.200 

7.5351 

0.3103 

7.6245 

720. 

1333.440 

4.96  39 

0.1025 

-6.6188 

WIND, 

1. 

HOUR 

1.852 

20 

5.8491 

7 

4.5236 

-58.1521 

Interpolation  between 

2. 

3.704 

15.80  31 

25. 0616 

-22.8319 

Snapshots  4  &  6 

3. 

5.556 

31.6745 

31.1321 

-4.0379 

4. 

7.408 

45.2941 

44  .  7367 

3.0175 

5. 

9.260 

52.8056 

52. 1772 

6.3979 

6. 

11.112 

55.67  64 

54. 9776 

8.3532 

7. 

12.964 

55.5360 

55.7612 

9.8517 

8. 

14.816 

55.2909 

55.4289 

11.1628 

9. 

16.668 

55.4900 

54. 5386 

12.3841 

10. 

18.520 

54.2968 

53.2621 

13.5262 

15. 

27.780 

43.6838 

47.4276 

18.2287 

20. 

37.040 

45.12  81 

44.1740 

25.3132 

.  25. 

46.300 

42.5121 

41. 9100 

32.3032 

30. 

55.560 

43.5284 

40  o 1261 

37.4110 

40. 

74.080 

35.15  86 

35. 9069 

42.6494 

50. 

92.600 

31.4613 

31. 2270 

43.5405 

60. 

111.120 

27.2626 

27. 0013 

42.3935 

75. 

138.900 

22.2634 

21. 9178 

39.4404 

90. 

166.680 

13.  58  53 

18. 1153 

35.9769 

120. 

222.240 

13.8112 

12. 9681 

27.9855 

180. 

333.360 

3.1801 

7. 2034 

16.8972 

240. 

444.480 

7.92  65 

4. 0666 

12.5087 

360. 

666.720 

7.62  58 

1. 5307 

12.1834 

480. 

888.960 

7.57 13 

0.7078 

12.4845 

600. 

1111.200 

7.5267 

0. 3684 

8.9790 

720. 

1333.440 

4.95  87 

0. 1301 

-2.1435 

Figure  Cl  4.  (Sheet  4  of  12) 


Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gilbert 


C35 


- Wind  Speed -  Inflow  Anglo. 


Radius 

Radius 

Scalar  Avg. 

Vector  Avg. 

( +  -  in,  -  ■  a 

(n.m.) 

(km) 

(m/sac) 

(m/sec) 

(deg) 

WIND 

•  HOUR 

20 

8 

1. 

1.852 

7.0426 

5. 0716 

-54.7408 

2. 

3.704 

17.1757 

16. 5115 

-18.2037 

3. 

5.556 

32.48  57 

31. 9070 

-1.9445 

4. 

7.408 

44.2867 

43. 6128 

3.8353 

5. 

9.260 

53.55  59 

49. 7581 

6.2712 

6. 

11.112 

53.2640 

52. 4012 

7.5746 

7. 

12.964 

54.2964 

53. 3844 

8.7047 

8. 

1*. 816 

54.46  36 

53. 4694 

9.9306 

9. 

16.668 

54.1409 

53. 0254 

11.2704 

10. 

18.520 

53.46  84 

52. 2117 

12.6752 

15. 

27.780 

53.3944 

48. 7947 

20.3394 

20. 

37.040 

43.4448 

47.  5966 

31.2649 

25. 

46.300 

45.0761 

45. 6723 

39.2970 

30. 

55.560 

43.35  28 

43. 0839 

43.4863 

40. 

74.080 

35.8795 

36  .  67  2  3 

45.6317 

50. 

92.600 

31.0389 

30. 8200 

44.6116 

60. 

111.120 

25.3489 

26. 0853 

42.2429 

75. 

138.900 

21.2724 

20. 9049 

33.3987 

90. 

166.680 

17.7416 

17. 2313 

34.5203 

120. 

222.240 

13.2931 

22.  3823 

26.4577 

180. 

333.360 

9. 0301 

6. 9736 

16.3809 

240. 

44*.48G 

7.9248 

4. 0382 

12.9052 

360. 

666.720 

7.6209 

1. 6167 

12.8425 

480. 

888.960 

7.5641 

0. 7856 

13.0809 

600. 

1111.200 

7.5185 

0.4267 

9.9572 

720. 

1333,440 

4.9536 

0. 1582 

0.7499 

MIND, 

HOUR 

20 

9 

1. 

1.852 

5.99  81 

5.0795 

-54.8159 

2. 

3.704 

17.1768 

16.5351 

-18.2066 

3. 

5.556 

32.50  58 

31. 9429 

-1.9150 

4. 

7.408 

44.2801 

43. 6174 

3.8725 

5. 

9.260 

53.47  46 

49. 6803 

6.2548 

6. 

11.112 

53.1378 

52. 2742 

7.4287 

7. 

12.964 

54.1663 

53. 2625 

8.4029 

8. 

14.816 

54.3547 

53. 3862 

9.4955 

9. 

16.668 

54.0787 

53. 0050 

10.7557 

10. 

18.520 

53.4631 

52. 2558 

12.1347 

15. 

27.780 

53.72 77 

49. 1375 

19.9558 

20. 

37.040 

49.04  89 

48. 2190 

31.2943 

25. 

46.300 

45.7582 

46. 3703 

39.4880 

30. 

55.560 

44.0145 

43. 7565 

43.6872 

40. 

74.080 

37.40  74 

37. 2072 

45.7600 

50. 

92.600 

31.4470 

31. 2338 

44.6753 

60. 

111.120 

25.  66  54 

26. 4072 

42.2520 

75. 

138.900 

21.4959 

21. 1338 

38.3407 

90. 

166.680 

27.9024 

17. 3978 

34.4100 

120. 

222.240 

13.3829 

12. 4762 

26.2241 

180. 

333.360 

9.06  83 

7. 0029 

15.9467 

240. 

444.480 

7.95  28 

4. 0419 

12.4316 

360. 

666.720 

7.6355 

1. 6109 

12.3296 

480. 

888.960 

7.5744 

0.7805 

12.4548 

600. 

1111.200 

7.52  58 

0.4218 

9.4619 

720. 

1333.440 

4.95  65 

0. 1548 

0.8799 

Snapshot  6 


Interpolation  botwaan 
Snapshots  5  A  6 


C36 


Appendix  C  Sample  Application.  12  Snapshots  of  Hurricane  GHbert 


Redkis 

ftsdkts 

- Wind  Speed  — - 

Scalar  Avg-  Vector  Avg. 

Inflow  Angk 

(  +  -  in,  -  - 

(rum.) 

(km) 

(m/sec) 

(m/sec) 

(dag) 

WIND, 

1. 

HOUR 

1.852 

20 

S.9S  37 

10 

5. 0869 

-54.8961 

Cm 

3.704 

17.1782 

16. 55BS 

-18.2091 

3. 

5.556 

32.5267 

31.  9786 

-1.8853 

4. 

7.408 

44.27  50 

43. 6215 

3.9100 

5. 

9.260 

53.3961 

49.6013 

6.2389 

6. 

11.112 

53.0167 

52. 1447 

7.2816 

7. 

12.964 

54.0441 

53. 1380 

8.0985 

8. 

14.816 

54.2563 

53. 3016 

9.0584 

9. 

16.668 

54.0279 

52. 9847 

10.2410 

10. 

18.520 

53.4681 

52.  3011 

11.5963 

15. 

27.780 

51.06  40 

49. 4815 

19.5878 

20. 

37.040 

49.65  22 

48.  8398 

31.3240 

25. 

46.300 

47.4393 

47 . 0666 

39.6742 

30. 

55.560 

44.67  52 

44.  4276 

43.8827 

40. 

74.080 

37.9347 

37. 7411 

45.8854 

50. 

92.600 

31.85  48 

31. 6*70 

4*. 7490 

60. 

111.120 

26.9819 

26.7288 

42.2611 

75. 

138.900 

21.71  96 

21.3627 

38.2847 

90. 

166.680 

13.06  38 

17. 5646 

34.3006 

120. 

222.240 

13.47  30 

12. 5704 

25.9949 

180. 

333.360 

9.1069 

7.0325 

15.5172 

240. 

444.480 

7.9816 

4. 0458 

11.9632 

360. 

666.720 

7.6511 

1. 6052 

11.8140 

480. 

888.960 

7.5856 

0.7755 

11.8209 

600. 

1111.200 

7.5338 

0.4170 

8.9545 

720. 

1333.440 

4.9598 

0. 1515 

1.0131 

WIND, 

1. 

HOUR 

1.852 

20 

3.7207 

11 

7.8357 

-13.2877 

2. 

3.704 

16.9519 

16. 5147 

-6.27B4 

3. 

5.556 

26.1567 

25  .  7  5  71 

0.5063 

4. 

7.408 

32.99  96 

32. 5674 

3.7029 

5. 

9.260 

36.6307 

36. 1683 

5.0808 

6. 

11.112 

33.31  53 

37. 6633 

5.7924 

7. 

12.964 

39.0787 

38. 6623 

6.5059 

8. 

14.816 

39.26  26 

38. 8826 

7.4000 

9. 

16.668 

39.96  61 

38. 7132 

8.4502 

10. 

16.520 

33.5591 

38  .  22  41 

9. 5350 

15. 

27.780 

35.5052 

36. 172* 

13.9470 

20. 

37.040 

35.8351 

36  .  5  3  06 

20.8314 

25. 

46.300 

33.2610 

37.9182 

28.2281 

30. 

55.560 

39.2255 

3B.8999 

34.0331 

40. 

74.080 

35.82  23 

36.  6041 

40.1687 

50. 

92.600 

32.2816 

32. 0769 

41.0231 

60. 

111.120 

23.2033 

27.9708 

39.6747 

75. 

138.900 

23.45  24 

23. 1454 

36.7160 

90. 

166.680 

23.0079 

19.6021 

33.6159 

12C). 

222.240 

15.5677 

14. 8986 

27.1520 

180. 

333.360 

11.10  93 

9.  7741 

18.6237 

240. 

444.480 

3.9368 

6. 8459 

15.5806 

360. 

666.720 

7.86  36 

3. 9343 

14.4370 

480. 

888.960 

7.6628 

2.5745 

14.8996 

600. 

1111.200 

7.5550 

1.8281 

14.1238 

720. 

1333.440 

4.95  76 

0.9003 

12.5738 

Figure  Cl  4.  (Sheet  6  of  12) 


Snapshot  6 


Interpolation  botwson 
Snapshots  6  Ml  7 


Appendx  C  Sample  Application,  12  Snapshots  of  Hurricane  Gilbert 


037 


-  Inflow  Anglo. 


Radius 

Radius 

Scalar  Aug. 

Vector  Avg. 

(+  ■  In,  •  ■ 

(tun.) 

flun) 

(m/sec) 

(m/sac) 

(dag) 

WIND 

.HOUR 

20 

12 

1. 

1.852 

12.3410 

11. 3551 

3.7061 

2. 

3.704 

17.37  32 

17. 0704 

5.5010 

3. 

5.556 

19.82  49 

19. 5672 

4.8894 

4. 

7.408 

21.3645 

21.  1058 

4.0851 

5. 

9.260 

22.32  64 

22. 0409 

3.5390 

6. 

11.112 

22.9808 

22. 6516 

3.4733 

7. 

12.964 

23.42  77 

23. 0439 

3. 8438 

8. 

14.816 

*3. 66  8 6 

23. 2223 

4.4889 

9. 

16.668 

23.7225 

23. 2100 

5.1796 

10. 

16.520 

23.58  59 

23. 0055 

5.6476 

15. 

27.780 

23.2624 

22. 4453 

1.8611 

20. 

37.040 

2S.1342 

25. 5009 

0.3562 

25. 

46.300 

31.03  52 

30. 4071 

10.2215 

30. 

55.560 

35.07  99 

34. 4948 

21.2548 

40. 

74.080 

35.06  79 

35. 7961 

34.1508 

50. 

92.600 

32.8376 

32. 6266 

37.4436 

60. 

111.120 

23.48  95 

29.  2604 

37.3449 

75. 

138.900 

25.21 07 

24. 9394 

35.4033 

90. 

166.680 

21.98 17 

21. 6405 

33.0628 

120. 

222.240 

17.70  00 

17. 1812 

28.0628 

180. 

333.360 

13.34  70 

12. 4114 

21.4032 

240. 

444.480 

13.91 23 

9. 5324 

17.4889 

360. 

666.720 

3.6205 

6. 1845 

15.2893 

480. 

888.960 

7.94  91 

4. 2965 

15.4385 

600. 

1111.200 

7.6724 

3.1780 

14.7235 

720. 

1333.440 

5. 00  44 

1.6352 

13.5201 

WIND, 

HOUR 

20 

13 

1. 

1.852 

7.7009 

6. 5613 

-2.0170 

2. 

3.704 

11.4624 

10. 8*61 

-1.3424 

3. 

5.556 

14.6535 

14.2437 

-0.4667 

4. 

7.408 

17.3708 

17. 0409 

0.1374 

5. 

9.260 

19.5312 

19.  2268 

0.5299 

6. 

11.112 

21.1956 

20. 8754 

0.9853 

7. 

12.964 

22.4706 

22. 1060 

1.7251 

8. 

14.816 

23.38  62 

22. 9570 

2.6868 

9. 

16.668 

23.9666 

23. 4606 

3.6873 

10. 

18.520 

24.2226 

23.  6335 

4.4648 

15. 

27.780 

24.31 32 

23. 9490 

1.5678 

20. 

37.040 

23.1526 

27. 4695 

0.8347 

25. 

46.300 

33.2867 

32. 4993 

10.3299 

30. 

55.560 

37.88  70 

37.  1112 

21.6899 

40. 

74.080 

43.1685 

39. 8787 

37.0295 

50. 

92.600 

37.22  99 

37. 0326 

41.6416 

60. 

111.120 

33.64  95 

33. 4597 

42.3094 

75. 

138.900 

23.6322 

28. 4091 

40.7412 

90. 

166.680 

24.6445 

2*. 3656 

38.3392 

120. 

222.240 

19.16  97 

18. 7279 

32. 5782 

180. 

333.360 

13.62  40 

12.7291 

23. 5062 

240. 

444.480 

13.6611 

9. 2116 

17.66B3 

360. 

666.720 

8.2944 

5. 3560 

14.5539 

480. 

888.960 

7.78  99 

3. 4564 

14.7964 

600. 

1111.200 

7.5948 

2.4277 

14.0900 

720. 

1333.440 

4.97  36 

1. 1949 

12.8509 

Snapshot  7 


Interpolation  between 
Snapshots  7  &  8 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gibert 


- Wind  Spaed - 

Inflow  Angk 

Radius 

Radiua 

Scalar  Avg. 

Vector  Avg. 

(♦  -  in,  -  - 

liun.) 

(km) 

(m/sec) 

(m/aac) 

(dag) 

WIND, 

HOUR 

20 

14 

1. 

1.852 

5.7881 

1.0245 

-54.9114 

2. 

3.704 

5.4541 

3.  9720 

-24.5178 

3. 

5.556 

9.4686 

8.6090 

-11.9362 

4. 

7.408 

13.4829 

13. 0106 

-6.3941 

S. 

9.260 

15.7968 

16. 4585 

-3.6516 

6. 

11.112 

19.4415 

19.1273 

-2.0103 

7. 

12.964 

21.5386 

21. 1914 

-0.5934 

S. 

14.816 

23.12  55 

22.7110 

0.8438 

9. 

16.668 

24.2274 

23. 7256 

2.2342 

10. 

18.520 

24.86  98 

24. 2708 

3.3539 

15. 

27.780 

25.3710 

25. 4490 

1.3220 

20. 

37.040 

33.1687 

29. 4218 

1.2605 

25. 

46.300 

35.5447 

34. 5681 

10.4416 

30. 

55.560 

43.6915 

39. 6964 

22.0873 

40. 

74.080 

44.26  57 

43. 9524 

39.3821 

50. 

92.600 

41.6690 

*1. 4767 

44.9388 

60. 

111.120 

37.  90  86 

37. 7362 

46.1393 

75. 

138.900 

32.1717 

31. 9795 

44.8648 

90. 

166.680 

27.4255 

27. 1891 

42.4840 

120. 

222.240 

23.72  41 

20. 3404 

36.3204 

180. 

333.360 

13.90  54 

13. 0482 

25.5059 

240. 

44h.480 

13.4115 

8. 8880 

17.8943 

360. 

666. 720 

3.0717 

4. 5179 

13.6535 

480. 

888.960 

7.6928 

2.6034 

13.7065 

600. 

1111.200 

7.S5  79 

1.6662 

12.8807 

720. 

1333.440 

4.9537 

0. 7413 

11.2863 

HIND, 

HOUR  . 

20 

15 

1. 

1.852 

5.27  86 

3.6952 

-21.4673 

2. 

3.704 

9.6598 

8. 8585 

-8.8558 

3. 

5.556 

13.9015 

13.4537 

-3.6465 

4. 

7.408 

17.2845 

16.9501 

-1.3335 

5. 

9.260 

19.90  31 

19.5926 

-0.1562 

6. 

11.112 

21.8548 

21. 5198 

0.7374 

7. 

12.964 

23.2916 

22. 9027 

1.7017 

8. 

14.816 

24.2837 

23  .  82  0  7 

2.7128 

9. 

16.668 

24.8923 

2*. 3430 

3.6006 

10. 

18.520 

23.1624 

24. 5209 

4.1324 

15. 

27.780 

25.32  68 

25.4590 

0.4130 

20. 

37.040 

33.7382 

30. 0349 

2.6953 

25. 

46.300 

35.3498 

35. 5371 

14.4592 

30. 

55.560 

43.5023 

39.6200 

25.8454 

40. 

74.080 

41.3230 

41. 2674 

39.3602 

50. 

92.600 

37.98  30 

37.7975 

43.0899 

60. 

111.120 

34.0207 

33.8370 

43.3371 

75. 

138.900 

23.6710 

26. 4499 

41.4519 

90. 

166.680 

24.47  93 

24.1978 

38.8298 

120. 

222.240 

13.7628 

18. 3017 

32.6069 

180. 

333.360 

13.0273 

12. 0419 

22.6556 

240. 

444.480 

13.0310 

8.3815 

16.5407 

360. 

666.720 

S.0627 

4. 5075 

13.6395 

480. 

888.960 

7.6978 

2.7350 

13.7344 

60  0  . 

1111.200 

7.55  89 

1.8231 

12.9673 

720. 

1333.440 

4.95  34 

0.8437 

11.6417 

Snapshot  • 


Interpolation  batwoon 
Snapahota  8  ft  9 


Appendix  C  Semple  Application,  12  Snapehota  of  Hurricane  Gifcert 
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- Wind 

Spaed  — - 

Inflow  Angh 

Radius 

Radius 

Scalar  Avg, 

Vector  Avg. 

(  +  -  In,  -  -  i 

(run.) 

(km) 

(m/sac) 

(m/sac) 

(dag) 

WIND. 

HOUR 

20 

16 

1. 

1.852 

7.5921 

6. 3574 

-15.8577 

2. 

3.704 

13.8447 

13.4131 

-2.9746 

3. 

5.556 

13.2339 

17. 9469 

0.7354 

4. 

7.408 

21.1606 

20.8  9  2  5 

1.7934 

5. 

9.260 

23.0540 

22.7  5  3  8 

2.2868 

6e 

11.112 

24.2866 

23.9272 

2.8904 

7. 

12.964 

23.0617 

24. 6285 

3.6555 

8. 

14.816 

25.45  96 

24. 9436 

4.4044 

9. 

16. 6o8 

25.5713 

24. 9689 

4.8958 

10. 

18.520 

25.4633 

24. 7739 

4.8956 

15. 

27.780 

25.2953 

25.4727 

-0.4949 

20. 

37.040 

31.3771 

30. 6556 

4.0804 

25. 

46.300 

37.3395 

36.6335 

18.2633 

30. 

55.560 

43.5397 

40. 0771 

29.5947 

40. 

74.080 

33.7619 
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Appendix  C  Sample  Application,  12  Snapshots  of  Hurricane  Gribert 
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- Wind 

flpud  — — — 

Inflow  Anglo. 

Radius 

Radius 

8oaiar  Avg. 

V actor  Avg. 

(♦  m  in.  -  m  out) 

(run.) 

(km) 

(m/aao) 

(m/aao) 

(dag) 

WIND, 

HOUR 

20 

22 

I. 

1.852 

5.82  55 

0. 2000 

-87.6107 

2. 

3.70* 

5.8812 

0.5546 

-70.1721 

3. 

5.556 

5.0218 

1.8204 

-37.1198 

4. 

7.408 

5.7394 

4.6981 

-21.8820 

5. 

9.260 

9.6611 

8.8197 

-14.1624! 

6. 

11.112 

13.5649 

13.0636 

-9.2704 

7. 

12.964 

17.2198 

16.8280 

-6.1336 

8. 

14.816 

21.57  20 

20. 1906 

-3.9692 

9. 

16.668 

23.48  62 

23. 0649 

-2.2587 

10. 

18.520 

25.37  45 

25.3865 

-0.7408 

15. 

27.780 

32.5430 

31.8957 

3.2647 

20. 

37.040 

35.4582 

34. 8746 

7.5786 

25. 

46.300 

35.5038 

35  .  8823 

13.0161 

30. 

55.560 

37.1546 

36.4978 

IB. 5939 

40. 

74.080 

35.71  30 

36. 3500 

29.7262 

50. 

92.600 

34.3014 

34.0590 

35.2759 

60. 

111.120 

31.4108 

31.1831 

36.9782 

75. 

138.900 

27.08  41 

26. 8226 

36.4202 

90. 

166.680 

23.4467 

23. 1226 

34.5579 

120. 

222.240 

13.2176 

17.7032 

29.1176 

180. 

333.360 

12.62  31 

11.5294 

19.5335 

240. 

444.480 

9.6441 

7. 7750 

13.5450 

360. 

666.720 

7.9456 

3. 8433 

11.0326 

480. 

888.960 

7.65  50 

2. 1516 

10.7810 

600. 

1111.200 

7.5351 

1. 3193 

9.9973 

720. 

1333.440 

4.94  59 

0. 5603 

8.7143 

FORTRAN  STOP 


Figure  Cl  4.  (Sheet  12  of  12) 


Snapshot  12 


Appendx  C  Sample  Appficaion.  12  Snapahofc  of  Hurricane  Gdbert 


C43 


Appendix  D 
Sample  Application  of 
Upgraded  CE  Model  to 
Simulation  of  36-Hr  Period 
of  Hurricane  Gilbert  in  the 
Gulf  of  Mexico 


This  appendix  provides  information  related  to  a  36-hr  simulation  of  Hurri¬ 
cane  Gilbert  in  the  Gulf  of  Mexico.  The  simulation  was  performed  by  OWI  as 
an  additional  test  with  the  upgraded  CE  model  The  simulation  time  period 
begins  1200  UTC  (Universal  Time  Coordinate,  formerly  known  as  Greenwich  ' 
Mean  Time)  IS  September  1988.  The  four  snapshots  used  in  the  simulation 
all  took  advantage  of  the  double  exponential  form  for  pressure  profile  specifi¬ 
cation.  Input  file  information  on  the  snapshots  and  storm  track  specification  is 
provided  in  this  appendix. 

The  appendix  also  includes  {dots  of  Hurricane  Gilbert  wind  fields.  Model 
wind  fields  at  19-m  elevation  are  given  at  6-hr  intervals  throughout  die  simula¬ 
tion  period.  Wind  speed  and  direction  is  represented  with  the  conventional 
weather  map  “wind  barb”  notation.  The  shaft  of  each  wind  arrow  indicates  the 
direction  and  die  barbs  or  “feathers”  indicate  wind  speed.  A  half-barb  denotes 
5  knots,  full  barb  denotes  10  knots,  and  solid  flag  denotes  SO  knots. 
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Figure  D1 .  Snapshot  and  storm  track  specification,  Hurricane  Gilbert 
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Figure  D2.  Modelled  surface  wind  field  in  Hurricane  Gibert  (Sheet  1  of  7) 
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Figure  D2.  (Sheet  2  of  7) 
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13.  (Conducted). 


The  upgraded  model  is  demonstrated  with  historical  hurricanes.  The  five-nest  and  seven-nest  models  are  applied 
to  Hurricane  Camille.  The  fully  upgraded  model,  with  seven  nests  and  general  pressure  specification,  is  applied  to 
Hurricane  Gilbert.  This  hurricane  was  chosen  because  it  is  well-documented  by  Black  and  Willoughby  (1992)  and  it 
evolved  into  some  nontraditional  storm  structures.  The  upgraded  model  was  more  effective  than  the  standard  CE 
model  in  simulating  the  storm. 


